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1 Introduction 



This paper is the fourth in a series devoted to numerical evaluation of the multi-loop, multi-leg 
Feynman diagrams that appear in any renormalizable quantum field theory. In |1| (hereafter I) the 
general strategy has been designed and in ^ (hereafter II) a complete list of results has been derived for 
two- loop functions with two external legs, including their infrared divergent on-shell derivatives. Results 
for one- loop multi-leg diagrams have been shown in |5] and additional material can be found in |^. 

In a new series of papers we will present a complete set of results for the class of fully general three- 
point, two- loop diagrams. There has been an extensive writing on the subject and known results in this 
area are given in however, to the best of our knowledge, we are still missing some general attempt to 
systematize the field (with the noticeable exception of [6 ). 

Due to the complexity of the problem, six different topologies with up to six propagators, see Fig. [3 
we have found natural to split our presentation into three parts: in this paper we will analyze general, 
scalar, massive configurations, while the treatment of the corresponding tensor integrals ,7, and infrared 
divergent two-loop vertices |H] will be discussed in two forthcoming papers (for recent numerical tests of 
existing analytical results see jH] and JO]). Furthermore, there are parallel developments in our project: 
a FORM [T^ code is being created for the generation of all one- and two-loop diagrams (so far only in 
the context of the standard model [T^), and all the numerical algorithms are being assembled into a 
FORTRAN code ^]. Well-known packages for diagram handling are in |14j . 

Any Feynman diagram is built, starting from the rules fixed by the Lagrangian of the theory under 
study, using two ingredients: vertices and propagators. The latter are represented in momentum space 
by 

^(^) (1) 

_|_ j^Z — i() 

where 5 — > 0+ and N{p) is a structure depending on spin. For a selected number of external legs and a 
given order in perturbation theory, the diagrams are constructed from propagators and vertices according 
to the allowed topologies, and are then endowed with signs and combinatorial factors [El- 

Despite the dramatic recent progress in the field of analytical evaluation of Feynman integrals, 
Mellin-Barnes techniques JHI) shifting dimensions differential ^Hl and difference ^HI equations, 
non-recursive solutions of recurrence equations f^ni, just to name a few, we firmly believe that sooner or 
later we will achieve the structural limit of the field and that a general solution can only be numerical. 

Our general approach toward the numerical evaluation of an arbitrary Feynman diagram G, a typical 
example of a multi-scale problem, is to use a Feynman parameter representation and to obtain - diagram- 
by-diagram - some integral representation of the following form: 

G = ^ [ dxQix), (2) 

Js 

where x is a vector of Feynman parameters, S is some simplex, Q is an integrable function (in the limit 
6^0) and Bq is a function of masses and external momenta whose zeros correspond to true singularities 
of G, if any. The Bernstein-Tkachov (hereafter BT) functional relations |2j are one example but, in 
this work, we will consider other realizations of Eq.@. One example will be enough to characterize the 
criterion of smoothness. Although the integral 



/= / dx 



V{x)-i6^ =J dx{hx^ + 2kx + l-i5j (3) 

is well defined through the i 5 prescription, it is not convenient to attempt a direct numerical integration 
because poles may be dangerously close to the real axis. However, after some algebraic manipulation we 
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can write Vs = V — i5, 



I = 



hl-k^ 



h 




(4) 



which is much easier to evaluate and where hi = k'^ represents a pinch singularity |^ for I. One 
drawback of Eq.Q is that the behavior around the singularity is overestimated. Therefore, whenever in 
the presence of a true singularity of G, we took particular care in investigating the integral representation 
with the correct behavior since, as well-known and as shown in our specific example, most if not all of 
the methods currently employed tend to overestimate the singular behavior. 

Smoothness requires that the kernel in Eq.((2)) and its first A'^ derivatives be continuous functions and, 
ideally, should be as large as possible. However, in most of the cases we will be satisfied with absolute 
convergence, e.g. logarithmic singularities of the kernel. This is particularly true around the zeros of Bq 
where the large number of terms required by obtaining continuous derivatives of higher order leads to 
large numerical cancellations. 

There are other alternative approaches to the numerical evaluation of multi-loop Feynman diagrams 
where differential equations are written and solved but, to our knowledge, they have been applied - so 
far - only to two-point functions or to QED/QCD examples |23j . 

In any realistic calculation tensor integrals will appear and we distinguish between two different 
applications: before assembling all terms into the 5-matrix element associated with some physical process 
we must verify the complete set of Ward-Slavnov- Taylor identities (hereafter WST) and this requires 
reduction, often called scalarization, of the corresponding expressions. While at one-loop level we only 
need reduction to standard scalar functions Ao,Bo, ■ ■ ■ j24j . here the number of scalar expressions which 
are required is larger because of irreducible numerators, i.e. we have less propagators in a diagrams than 
scalar products, and the adjective "scalar" should be used only when we enlarge the class of functions to 
arbitrary space-time dimension and arbitrary powers for the propagators. However, in [2| we will be able 
to show that this is only a semantic question and we will present our version for a scheme of scalarization. 
After checking that the relevant WST identities are satisfied, we plan to organize our calculation for 
physical observables according to gauge-parameter independent blocks Bi, each of which will be mapped 
into one integral of the form 



where P are polynomials in the Feynman parameters. Of course, checking unrenormalized WST identities 
is only a test on having correctly generated and evaluated the diagrams and, in the end, we will have to 
control also renormalized WST identities. 

We mention here that, in our approach, numerical evaluation of infrared divergent vertices is always 
understood as numerical evaluation of the residues of the infrared (or even collinear) poles at space- 
time dimension n = 4 and of the corresponding finite parts. In any realistic calculation infrared (and, 
eventually, collinear poles) will be combined with real emission of both photons and gluons according to 
some scheme, e.g. dipole- formalism [251 • strategy, to be fully described in |H1, will always start with 
techniques which aim to extract the infrared/collinear divergent part of two-loop vertices, with a finite 
remainder. To summarize we either express the diagram through hypergeometric functions in arbitrary 
space-time dimensions that are subsequently expanded around n = 4 (cf. with Sect 7.4 of |2j) or we use 
sector decomposition (a technique introduced in [201 )• both cases we are able to obtain results valid 
in all regions, not only in the unphysical one. 

Finally we stress that, whenever possible, we have been seeking for at least two independent algorithms 
per diagram in order to allow for an internal cross-check of our results. 




G 



(5) 
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The outline of the paper will be as follows: in Section |21 we describe our conventions for a general 
definition of two-loop diagrams. Special cases and special topologies are discussed in Section |31 The class 
of G^^^ diagrams is analyzed in Section |3J All two-loop vertex-like topologies and the related derivatives 
of self-energy diagrams needed in the wave function renormalization factors are described and evaluated 
in Section \^ - IIUI Numerical results are shown in Section 1111 Several technical details, ranging from 
solutions of Landau equations |27j to recurrent integrals encountered in the text are analyzed in several 
appendices. 



2 General definitions and conventions 

In this section we present basic definitions and properties of Feynman diagrams. An arbitrary two-loop 
scalar diagram can be cast in the following form: 

(a; mi, . . . , mo,; r/i|7; rria+i, . . . , ma+^; 7/i2|/3; rria+^+i, . . . , ma+^+/3; 7^2) = 
'"^ J d-q.d-q^llik^+m^)-' II {k]+m])-' H {kf + mfr\ (6) 

i=l j=a+l i=a+-y+l 



where n = 4 — e, n is the space-time dimension, and a, f3 and 7 give the number of lines in the qi,q2 and 
Qi — Q2 loops, respectively. Furthermore we have 

ki = qi + EjLi ^ij Pj^ i = l,...,a 

h = qi-q2 + Y.j=i VijPj, i = a + 1, . . . , a + 7 

ki = q2 + J2j=i VijPj, i = a 7 + 1, . . . , a 7 + /?, 

where N is the number of vertices, r]"" = ±1, or 0, {p} is the set of external momenta and fi is the 
arbitrary unit of mass. Diagrams that can always be reduced to combinations of other diagrams with less 
internal lines will never receive a particular name. Otherwise we will denote a two-loop scalar diagram 
with G"^'^ where G = S,V, B etc. stands for two-, three-, four-point etc. Family of diagrams with the 
same number of internal lines will be denoted collectively by S^, etc. 

Next we recall few basic properties of two-loop diagrams PF. In any two-loop diagram there are 
three one-loop sub-diagrams, called 07, f3'y and a(3 sub-diagrams respectively. 

Definition 1: the aP'j diagram is overall ultraviolet (UV) convergent if a -|- /? -|- 7 > 4, logarithmically 
divergent ii a + P + j = A, linearly divergent if a + (3 + ^ = 3, and so on. 



Definition 2: the q/3 sub-diagram is convergent if a + P > 2, logarithmically divergent if a + P = 2, 

3 

2' 



linearly divergent if a -|- /? = |, and so on. 



If any of the one-loop sub-diagrams diverge, we will have counter-terms associated with them. Therefore, 
in addition to the aPj diagram we will also consider the subtraction diagrams of Fig. □ G^l etc. More 
specifically, an arbitrary vertex diagram (often referred to as a three-point function) is depicted in Fig. [21 
where, in our conventions, all external momenta are flowing inward and P = P1+P2 ^- The sign of the 
invariants is left unspecified and each external mass squared can represent not only a mass squared but 
one of the Mandelstam invariants, s > and t < 0. We also introduce scaled quantities 

2 2 
^*' = ^^' ^ = 1>---^' ^1=1^1' ■?' = 1'2, (7) 



^In our metric, space-like p implies = + p\> 0. Also, it is pi = ipo with po real for a physical four-momentum. 
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etc. 



Figure 1: The arbitrary two- loop diagram Gt^^ of Eq.® and one of the associated subtraction sub-diagrams. 



where m^, i = 1, . . . ,N are the internal masses and Mandelstam invariants 



pf = -SiMf, Sp,Si = ±l, 



(8) 



where and are positive. The diagrams are evaluated within dimensional regularization where the 
space-time dimensionality is n = 4 — e. Therefore, we define a two-loop MS factor 

Ac,v = T + lnvr-Mn— , (9) 

where 7 = 0.577216 • • • is the Euler's constant. We will also use the shorthand notation 

S{n) = (^) ' r (n + 6) , G{n, M') = (-^) ^ T {n + e) , (10) 

where T denotes the Euler's gamma-function. Note that in this paper we will use the van der Bij-Veltman 
parametrization [25] of two- loop integrals rather than the more familiar Cvitanovic-Kinoshita one j3Uj . 
Finally, to keep our results as compact as possible, we introduce the following notations where xq = yo = 1 




Figure 2: The three-point Green function. AU external momenta are flowing inward, P — pi + P2- 



dSn{{x}) f{xi, ■■■ ,Xn) = '[]_ 



dCn{{x}) f{xi, ■■■ ,Xn) 

dCS{{x}; {y}) /(xi, • • • , a^^^, yi, • • • , y^J 



i=l 



dxi /(xi, ■ ■ ■ , Xji), 



/ TT /(xi,---,Xn), 
■^0 i=l 



1 '«1 



"2 



u i=l j=l -^u 



Vj-i 



dyj f{xi,- 



(11) 
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When no confusion arises we will, eventually, omit the list of arguments which will be, otherwise, made 
compact by using the following notation: 

{f{{x} ; y,z) for i = 
/({x}; fori = l . (12) 
f{{x} ; z,u) for i = 2 

Furthermore the so-called '+'-distribution will be extensively used, e.g. 



lo X J Jo X 

etc. Gram's determinants associated with N vectors Pi, • • • ,P2 are always denoted by 

Gi...t, = -det{pi-pj). (14) 

Quite often we are interested in the condition G12 > which requires that pi and/or p2 and/or pi + P2 
be time-like. Therefore the condition is violated only if the momenta occurring in G12 are all space-like, 
indeed time-like momenta and G12 < cannot happen for pi^2 real. 

2.1 Alphameric classification of diagrams 

As we have seen any scalar two-loop diagram is identified by a capital letter (S, V etc) indicating the 
number of external legs and by a triplet of numbers (a, /3 and 7) giving the number of internal lines (in 
the qi,q2 and qi — 52 loops respectively). There is a compact way of representing this triplet: assume 
that 7 7^ 0, i.e. that we are dealing with non-factorisable diagrams, then we introduce the integer 

K = 7max amax {(3 - I) + a - I +7 (15) 

for each diagram. For G = ^ we have amax = 2 and 7max = 2. We can then associate a letter of the 
alphabet to each value of k, thus introducing the following correspondence between a/37 and k: 

121 E, 131 I, 141 M, 221 G, 231 K, 222 H. (16) 

This classification is extensively used throughout the paper. 




Figure 3: Examples of two-loop vertices that are topologically equivalent to self-energies. 



3 Special kinds of diagrams 

In this section we introduce special diagrams, topologically equivalent to either one-loop or two-loop 
configurations, which are naturally connected with two-loop vertices. First of all we have two- loop vertices 
with the topology of a self-energy (i.e., diagrams with only two vertices connected to the external legs) 
which are shown in Fig. |31 They are given by S^'^^{P), S^^^{P) and S'^^^(P), with P = pi+p2-, and have 
been evaluated in (2j. 
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3.1 One-loop counter-terms in two-loop diagrams 



Renormalization in quantum field theory is a order- by-order procedure, therefore we also have to 
consider those one-loop vertices with the insertion of a counter-term, which are associated with the two- 
loop vertices of Fig.[3 There are several subtraction diagrams, for instance three associated with one- loop 






a) V^'^ Pk^ b) y^^^ \ c) V^'^ 

Figure 4: Examples of one- loop vertices with counter-terms (gray circles). Permutations are understood, 
vertex counter-terms. Fig. 0fa), 

1 



V 



sub 



(g2 + ml) {{q + pi)"^ + mf) {{q + pi + pj] 



m 



2\ ' 



(17) 



where i = 1, . . . , 3, and three associated with one- loop self-energy counter-terms. Fig. Efb), 



sub 



i TT^ 



(g2 + ml) {{q + pi)"^ + mf)^ {{q +Pi+ pj)'^ + 



m 



2\ ' 



where i = 1, . . . , 3. We easily obtain 



sub 



dVi Co(l,l, l\pl,p'j,pl;mk,mi,mj), F/.^" = 6SiCo{l,2,l\pi,pj,pi.;mk,mi,mj) 



sub 



„2 „2 2. 



(18) 



(19) 



where 6Vi,5Si have been fixed, for instance in the MS-scheme, by a one-loop calculation and contain an 
ultraviolet pole. Co is the generalized, scalar, form factor with arbitrary powers in the propagators which 
must be evaluated up to O (e), contrary to the one- loop case. According to the findings of J all one-loop 
diagrams are evaluated, at any order in e, according to the BT-algorithm. Counter-terms associated with 
one- loop four-point vertices. Fig. ^c), will be included as well: 



sub 



v;jk 



SV.'^Bo{pi;mi,m2). 



(20) 



Note that subtractions associated with the diagrams of Fig. |21 have not been explicitly included in Fig. 0] 



3.2 Wave function renormalization 

To perform two-loop wave function renormalization we need the derivatives of two-point functions 
which are always a combination of scalar and tensor self-energies (those with powers of irreducible scalar 
products in the numerator) where one of the propagators containing the external momentum p has power 
—2, as depicted in Fig. [3 



7 




Figure 5: The two- loop topologies that contribute to wave- function renormalization (permutations are not included). 
The 'o' refers to propagators with non-canonical power —2. 

All the diagrams of Fig. El represent special cases of two-loop vertices with one external momentum 
set to zero. By referring to Fig. [7] we obtain the following set of equivalences: 

5III = y 121 = ; ms = m^) , S^^^ = F^^i = ; m4 = mg) , 

^131 ^ y 141 = ; m4 = ms) , 5f ^ ^ F^si ^ g . ^ ^ (21) 

where Sp is the self-energy type of diagram where the propagator containing momentum p has non- 
canonical power —2. The previous relations have been derived for the scalar case, but they hold for 
tensor integrals too. The infrared divergent configurations for on-shell derivatives of two-loop two-point 
functions have been previously considered in II while the infrared finite ones will be treated, in this paper, 
as a special case of vertices. 



4 The G^^^ class of diagrams 

G^^ ^ is a general class of two-loop diagrams with -|- 2 internal lines which are overall ultraviolet 
convergent for N > 2 and contain the 07 logarithmically divergent sub-diagram. We have 

/-I '^^'^ -1 
d''qirq2{ql+ml)-'[{q^-q2f + miy (fe + A:,,) V " , (22) 

i=0 

where the momenta ki are linear combinations of the external momenta Pj, ki = Pq + Pi + . . . + Pi. 
Often, we are going to use a special notation: pfj denotes the square of the difference of the four- 
momenta flowing through propagators i and j. Since three out of six ^/-topologies belong to this family, 
we found it appropriate do give a detailed description of its properties. In order to evaluate this specific 
class of diagrams we introduce Feynman parameters Zi for the propagators of the q2 loop and also a 
{zj}-dependent momentum and mass 

JV — 1 JV — 1 JV— 1 N — 1 

= Y[ dzi 5 {1 - Zi^) , Zl = ^^ Zi, P^l = ^^ ^ikifi, = Zi {kf + , (23) 

1=0 i=0 i=0 i=0 



8 



and obtain 



G 



IJVI 



-2g'^-^WdPI 



771 



1 4-e 



2 + e 



2-iV-e 



(24) 



where x and y parameters have been introduced to combine all propagators arising after the q2 integration. 
Additional quantities are as follows: 



X 



m 



2 (1 — x) m\ + X 



x{l — x) 



j d[P] = j dC2{x,y) J 



-e/2 



dC2ix,y) I dz^ X{l-X) ' y^-^l-yyr^. 



Y/2 



(25) 



(26) 



The G^^^ family has special properties and the corresponding diagrams are particularly easy to handle 
within the BT method. First of all, due to the presence of the factor + y^ we can apply for the 
first term of Ea.H24|) the following BT relation: 



K + y') x'-^-' = (i + ^^3^ y dy) X 



2-JV-e 



Henceforth, an integration by parts gives 



dyy-\l-yYl^ml + P^y^x'—^ 



^ ^ 2(A^-2 + e) 



(27) 



(28) 



It is easily seen that the term proportional to 4 — e in Eq. (|28|) cancels the equivalent term in Eq. (|24|) . 
Henceforth, using well-known properties of the Euler's T function we obtain the following compact result: 



G 



iNl 



-g{N - 2) dx dy I dz, x(l - x) ^'^ (1 - yf'^-' x'"^"^ 



(29) 



As expected, for N > 2 the diagram is overall ultraviolet convergent while the sub-divergence is hidden 
in the y-integration. We use the (5-function to carry through the integration over and introduce a new 
set of variables defined by 

Ui - Ui+l 



y = uq, 



Zi 



Uo 

obtaining < Ujv-i < • • • < ^^i < "Uq < 1- Furthermore: 



iNl 



g{N-2) / dx / dSN{uor--,U!,-i) x{l - x) {1 - uoY^^^' x^'"'' 



0, 



G 



where the quadratic form x bas been rewritten as x = u^'Hu + 2K,*u + C, 



m. 



i+3 



m: 



i+2 



There is a special case, Pq = = 0, where we redefine uq = y, and obtain 

G^^i = -g{N- 2) dx j dSNiy, ui,---, u^-i) [x{l - x)] (1 - y)^/2-i 



(30) 



(31) 



(32) 



(33) 



X = u^Hu + 2/C*ii + {ml - ml) {1 - y) + m|. 



(34) 
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where H etc are given by 
'Hij -- 

In this case the following relation holds: 



Pi ■ Pj, ICi — - {k} 



j+3 



m. 



1 



(y - 1)5, + 



1 



l,...,iV-l. 



(35) 



At+l 



(36) 



where C/* = —IC^TL ^ and 6 = rrig — /C*?^ ^/C. Alternatively we may introduce new variables v with 



u 



v-n-^JC, x = v^nv + {ml - ml) {I - y) + b, 



(37) 



where x is 11 ow an homogeneous form in v. 

For each topology there is a maximum value for N, denoted by A'^max; corresponding to every two- 
loop diagram G G^^^, i.e. Nyna.x{S) = 3, A^max(^) = 4, etc. They all correspond to an internal line 
with a one-loop self-energy insertion and we are always interested in the case where the masses in the 
two propagators adjacent to the insertion are equal, otherwise the diagram can be expressed through the 
difference of other (simpler) diagrams. 

There are alternative assignments of internal momenta: consider the case of V^^^, as given in Fig. [7| 
(d). With the momentum assignment of Fig. [J] we may have /cq = and 



ki=pi, k2=pi+P2, k3 = 0, or ki = 0, ^2=^1, k3=pi+p2. 
Accordingly, the momenta Pi,i = 1, . . . , 3 are 

Pi=pi, P2=P2, P3 = -P = -Pi - P2, or Pi = 0, P2=Pl, P3=P2, 



(38) 
(39) 



and H,K of Eq. (|35|) should be consistently evaluated. 

When N = A^^ax and the masses in the two propagators adjacent to the self-energy insertion are equal, 
we are in the special case where two propagators coincide. As a consequence, one of the u-integrations 
in Eq.(j2Sl) is trivial and can be shifted to the innermost position so that we can write 



dSK{y,ui, ■ ■ • , wjv-i) 



(40) 



in Eq.(|2Sl), with some important consequence to be discussed in the next section. 



4.1 Power-counting for G^^^ when 6 

The G^^^ diagrams are evaluated with repeated applications of Eg . (|36|) . When 6 — > 0, the diagram 
could show a sing ularity if < U^-i <■■■ <Ui<l, with U = -H'^ /C, see Ea.dH^. 

The BT procedure that we have described will overestimate the nature of the singularity^ and we 
have to find some general algorithm to derive the leading behavior of the diagram for 6^0. With A^ > 3 
we can write 



G 



iNl 



-Q{N-2) 







/ dx 


x{l — x) 


JO 



'IJVI 



7liVl 



dSNiy,uu---,u^-i) (l-yY^^'^ {u-U)'n{u-U)+K{l-y) + b 



2-N-e 



(41) 



Any other known procedure will also overestimate the nature of the singularity. 
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where the quantities U and b have been defined after Eq. H36|). while K = rn^ — m|. If the point P of 
coordinates {ui = Ui} is internal to the integration domain P, we can complete V to the unit hypercube 
since the integral will be regular in the complementary part [0, l]'^— P where the diagram can be computed 
with standard methods. Therefore we arrive at the following set of relations, where the integration is 
now over the hypercube: 



r>lNl 



J dCr.y'^'^-^ [iu-UYniu-U) + Ky + b 



2-N-e 



1 JV-1 2 



-If 

2 2 „i JV-1 



(n) 



dui 



fl'^"^ (u^Hu + Ky + bf- 



E ••• E / ^y[n(-l)"'«J""V du,]y^/^-'iQn,...n^.,+Ky + b) 

ni=l njv-i=l ° i=l ° 



\2-JV-e 



(42) 



We have introduced 



a 



(1) 



-Ui. 



JV-l 



o;^- Ui Uj . 



(43) 



In the following we suppress this string of indices in Q. Using the well-known Mellin-Barnes technique j31j . 
we obtain 



— / dsB{s,N -2 + e- s)p 

^ ^ J — ioo 



N-2+e-s 



I{sO, 



p = b-\ 



2 2 I N~l 1 

^(^)=E--- E / dy[lli-ir'a^^U du,]y^/''~'iQ + Ky) 

ni=l n]v_i=l ° i=l ° 



(44) 



This relation is valid for < Re s < — 2 + e {N > 3) and B denotes the Euler's beta-function. 

The /- integral appearing in Eo. 1)44(1 is computed via sector decomposition [^S]: it is the sum of N 
contributions, the first given by 

2 2 „1 iV-l y 

^°=E--- E / dy[lli-irat^ du;\y^/'-HQ + Ky)-^ 
1 1 J 1 J 



ni=l njv_i=l ^ i=l 

(iCjv(y,ni, • • • ,Mjv-i) y^"^""* 



{u^ -U^)n{u-U)y + K 



(45) 



where we performed the transformation Ui ^ y {ui — Ui)/a(^''''^ and we have been able to set e = 0. We 
obtain 

f\s) = -K-' j dCN-i2Fi{s,N-l-s;N-s;-^). (46) 

where = (ti* — C/*) 7i{u — U) and 2-^^i denotes the hypergeometric function. Since we are interested 
in the behavior of the diagram for \p\ — > oo we close the contour over the right-hand complex half-plane 
at infinity. The poles of the integrand are ats = A^ — 2 + A: with A: > integer; they are all simple but 
for k = \ which is double. From this we obtain 



Go 



const + 0{b), for 6^0. 



(47) 
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The remaining A'^ — 1 contributions for / of Eq. H44() have the following form: 



2 

,un-i) y ' 

n=l 



^ r 



(48) 



where we changed variables according to y — > Um y and Ui — > Um {ui - Ui)/af''^ for (i / m). We 
introduced: 



(49) 



We complete first the y-integration, obtaining as a result 



Jo 



(Q;^(")nm)-^Fi(s,-;l + -;-^ 



2 ' 



Using well-known properties of the hypergeometric function |31j we derive 



(50) 



2 r 



Jrr 



n=l 



S Jo A 

e/2-s n 



/ex T{s-e/2) 
\2) T{s) 



Q 



R{n) 
m 



N~2+e-2s 







(51) 



Note that Qm^"'^ is ifm-independent, allowing for Um integration. In particular, 
-1 

dumU^~'^~^ 2Fi{s , s; 1 + s; qnmUm) = B{l,N-l-s)3F2{s , s, N-l-s; s + l, N-1; g„„), (52) 







with qnm = Qrn^'^ /K and where 3-F2 is a generalized hypergeometric function. The general expression 
will be 



7liVl 



— (-1)" / dCN-2 / ds B{s, N - 2 - s) B{1,N -1-s 



fT'i^) pN-2-s 



V2/ r(s) ^ ^ A^-l + e-2s 



(53) 



where the functions ff^ are 



/r(s) = -E:-SF2(s, Af-l-s; s + l, iV-1; 



K 



QR(n) 



e/2-s 



■ (54) 



Closing once again the contour over the right-hand complex half-plane at infinity we have (for e = 0) 
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= 3: a single pole at s = 1, double poles at s = 2 + k, with k > for the first term in Ea. H53() : a 
double pole at s = 1 and single poles at s = 2 + A; (A; > 0) for the second one. 

> 3: a single pole at s = — 2 and double poles at s = N — I + k (k > 0) for the first term in 
Ea. (|5H|) : simple poles at s = (A^ — l)/2 and s = A^ — 2 + /c (/c > 0) for the second one. 

As a result, and taking into account the correct e dependence, we obtain 
Gm^ ~ ln6(ln^6) for divergent (finite) part, 



g^^^ ~ 5(3-JV)/2^^(3-Af)/2 ^^^^ ^ 3) divergent(finite) part. 



(55) 



for 6 — > 0. Ea.()55|) is derived under the condition that we have A^ + 2 different propagators in the diagram 
which is not the case if there is a self-energy insertion with equal masses in the two propagators adjacent 
to the insertion, e.g. for 5^^^ or V^'^^. When this happens we have to distinguish between Nj^ and Nj, 
i.e. the number of internal lines and the number of integration variables in Q^'^^ of Ea.()41|). For S^"^^ or 
y^^^ we have Nj = Nj^ — 1 since one integration is trivial and produces some extra factor in the integrand 
for Q^'^^ which, however, is irrelevant from the point of view of power counting. The most important 
consequence is that, in these cases, in Eq.(|53]) we will have 

1 



with a pole at s = A^/2 — 1 + e/2. Taking into account the ultraviolet pole of Eq. (|53j) we have a 
corresponding leading behavior of 

r(|) 61-^/2(1- I In 6). (57) 

Having discussed the simplest class of two-loop vertices in more detail, we now proceed to an exhaustive 
classification and description of all the vertices, from l/^^i y222^ 



5 The y^2i diagram 

The y^2i diagram of Fig. Ufa) is representable as 



TT V 



2e 



(58) 



[1] =ql+ rnl [2] = (gi - (72)' + ml [3] = {q2 - P2? 



m 



3' 



[4] = {q^-P)^+ml 



(59) 



Therefore it is a special case of G^^^, Eq. (|33j) . with ko = —p2 and ki = —p2 — pi- The discussion of the 
relative Landau equations is shown in Appendix IfI 



121 



5.1 Evaluation of V 

Following our discussion of G^^^ of Section |3] we are able to derive immediately the result for V^'^^. 
By performing an expansion in e we have: 



V, 



121 



g2 '^t/V 



121 

fin 



dCS {x ; y, z) 



1-2/ 



1-2 j dC2 lnxs(x,l,y) 
+ / dC2 lnx^(x,l,y)Li;^(x,y)-^-ic(2), (60) 
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where C is the Riemann's zeta- function, A^/v- is defined in Eq.@ and the subscripts DP, SP and fin 
denote the double and single ultraviolet pole (in dimensional regularization) and finite part (as n ^ 4), 
respectively, and where the auxiliary quantities are 



Li;Eix,y) = ln(l -y) -Inx - ln(l - x) - In Xglx, 1, y). 
For the complete list of different kinematic configurations we recall Eq.®: in general we have 

/ \ 22, 22,/ 2 2\ , r 2 2\ 



+ (S2 1^1 



fil)z + fil{l-y). 



(61) 

(62) 
(63) 

Note that Eg . (1601) defines V^'^^ everywhere. For the first integral in Eg . (16011 we may use the fact that, 
thanks to the i 6 prescription, both Xb(^' ^) ^^'^ Xsi^i 1' ^) have the same imaginary part and therefore 



We report the explicit result for the situation where Sp and si^2 are positive. Here we have 

/ N 22,22, /I 2 2\ /2,2 2\|/2 2\,2 

where we used the auxiliary quantities 

2 + 



/i?. = 1 + ^2 _ ^2 



lnXB(2;,y,2;) |+= X] ^"^2/ - 2/i 



(65) 



i=± 



where y±{x,z) are the roots of the equation Xsi-'^^y^ ^) ~ ^- ^^is way we obtain an expression which 
is more appropriate for numerical integration. 



dS2iy,z) 



lnXB(3;,y>2;) 



1-2/ 



E 



dz Li2 ( — - 
V 1 - 2/i 



(66) 



Li2 (z) denoting the standard di-logarithm 



111 



5.2 Wave function renormalization: S, 

The evaluation of the derivative Sp^^ follows immediately from Ea. l|6U|l with = 0, 
H3 = Hi, i.e. with 

Xe -^Xei^^y) =y^ - (1 + /^I - /^s) y + /^L 

where we assumed > 0. Therefore we have 

. O ^111 _ Am. ( . 1 



clll 
'-'p ; DP 



A 2 r,lll _2ri£^_L(A 



l-2 1n^i) 



1 and 
(67) 

(68) 



olll 
'-'p ; fin 



dC2y 



lnxg(3;, y) 
1-y 



+ 



dC2 ln^|Li;£;(2;,y) 



where L\-^e = lii(l — y) — Inx — ln(l — x) — In /ig. Finally we have 



rilll 
'-'p ; fin 



dC2y 



lnXB(a;,y) 



1-y 



■C(2), 



+ lnA.i(l-ln^i)----C(2). 



giving our result for the finite part of Sp^^. 



(69) 



(70) 
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6 The V^^^ diagram 

The V^^^ diagram of Fig. EJ^c) is representable as 

.'i'- = ,-/^„/<r.„^y^, ,71) 

with propagators 

[1] = ql + ml [2] = (gi - + ml [3] = + mf, 

[4] ^ {q2 + Pi? + ml [5] ^ (92 + Pf + ml (72) 

This topology represents a special case of G^'^^, Ea. (|5n|) . with = 0, ki = pi and k2 = Pi + P2- The 
corresponding Landau equations are discussed in Appendix [HI Two methods for the evaluation of V^^^ 
will be presented in the following sections. 

6.1 Evaluation of V^^^: method I 

Our first method for evaluating V^^^ uses the BT-method of Appendix^for increasing the powers in 
the integrand. Here we describe in more detail the derivation for time-like momenta. First we introduce 
the quadratic form 

/ \ 22, 22, /I 2 2\ ,^2,2 2\ 

+ (2^1-^45) ^2 + (M3-/^x)y + /^x) (73) 

with defined in Eq.ljTI) and fifj in Ea. (|BH) . Observing that Xii^^ lil/^-^) does not depend on x we will 
use XiiViZ) = Xi{x-,^,y, z). Furthermore, the coefficient hi is defined as 

bi = {^l + f^lf - ^35 (1 + - ^1) (^f + - mI) + ^! + Ai/ui, (74) 

where we have set Aj = A(l,i^^,t'|) and A denotes the familiar Kallen's lambda-function X{x,y,z) = 
x'^ + y'^ + z'^ — 2xy — 2xz — 2yz. BT co-factors fEa. (|269|) ) are 

Zo;i = -Xi, ^3;j = 0, Zi.j = {l-iyf-iy2){i^f-iJ.l) + 2{uf + f4-nl)u2 

Z2;i = -{l-vl- vl) {ul + fil-fil)-2 {vl - fil) vl = Zi-i - Zi+i-i. (75) 

An integration by parts is performed before the e-expansion and this leads to the following result: 

2 

Vg'' = {-- + ^uv)j^{2\i f dS2lnx,iy,z)+ j dCi ^Z^lnx,([l2/0],) + A,}. (76) 

The [x,y,z]i notation has been introduced in Eq.(|12|). Note that the expression multiplying the ultraviolet 
factor {AuY defined in Eq.®) has the form of a Co-function, the one-loop scalar vertex. For the ultraviolet 
finite part of the diagram we obtain the following expression: 



1/131 



M2 bi 



J dCS (x ; y, {z}) if + J dCS (x ; y, z) if + J dC2 I^ - ^ A,] . (77) 
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The various /-functions of Eq. (|77() are given by 



A, 



1-2/ 



i=0 



1 1 ^ 

-ff = A/ 1- lnx,(x,y,y,z) +lnx,(y,2) 1 - y, z) ]-^'^Zr^ 

1=0 



1-y 



(78) 



with an additional set of auxihary functions defined by 



Li-i{x,y,z) = C{x,y) - lnx^(y,^), L2-i{x,y) = C{x,y) - lnx^(l,y), 
L3;i{x,y) = C{x,y) - lnXj{y,y), Li.i{x,y) = C{x,y) - lnXj(y,0), 



(79) 



where C{x, y) = ln(l — y) — ln(x) — ln(l — x). Note that no modification of the algorithm is needed when 
the //-matrix becomes singular, i.e. for \j = 0. This method fails when bj, defined in Ea.(|74p. is zero. 
For a general configuration, Eq.®, we have 



Xj{x,y,Zi,Z2) = Si + 82^2 ^2 + (^p - SiV^ - 82^2)^1 ^2 " («! 1^1+^3- Hi) Zi 

+ {si vl - sp - /il + ^4) Z2 + (^3 - ^J'l)y + (so) 

Henceforth we seek for a second, alternative, algorithm which allows for internal cross-check. 

6.2 Evaluation of V^^^: method II 

Method I, described in the previous section, fails when we are around 6/ = or, equivalently around 



1/2 



2 4 = -ut 4) + + 4) + 4(1-4 + 4) ± [x{ut,4,4) x{i,4,4)^ ,(8i) 

which corresponds to the leading Landau singularity for all time-like momenta (other configurations follow 
with the appropriate change in signs). 

In we have given a complete discussion of the different options that one has in dealing with integrals 
where the related BT-factor is approaching zero. If we introduce Z = —Ti.~^ fC, with Tl, K, defined in 
Eq.JSSI), then 

X = A,,,{l-y) + {z-ZyH{z-Z) + h. (82) 

If 6/ = but the condition < Z2 < < 1 is not fulfilled, V^^^ is regular and we can perform a Taylor 
expansion around hi = 0. If, instead, the condition is satisfied, V^^^ is singular and we implement a 
Laurent expansion via Mellin-Barnes techniques following again pj. Also in jHj we have shown that for 
one- loop, multi-leg diagrams new integral representations can be constructed which encompass the need 
for expansion. 

In order to deal with y^^^ we start from the relation 



V 



131 



G{1) dCSix;y,{z}) x(l - x) (l-yy/'-^x 



-e/2 



Ne/2-1 



(83) 



where Q is defined in Eq. (|l()|) and we rewrite x ^ 

X = A^.j{l-y) + Bj{zi,Z2), Bj{zi,Z2) = z^nz + 21Ch + ml A^.j = ml-ml (84) 
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Now we introduce (always assuming time-like external momenta) 

Ax-i = nl - f^l, Bj = Xi{x,y,zi,Z2). 
The innermost integral in Ea. H83() will then become 

-l-e 



Ji 



l-Zl 



dyy 



e/2-l 



(85) 



(86) 



A similar result holds for all G^^^ diagrams. The integral appearing in Eg . (1861) is evaluated explicitly in 
Appendix [UJ The result is 



ln{l + Qx-i), Qx;i - 



A,/(l-^i) 



(87) 



where the Co-function is defined in Eg . (12821) and, in the limit e — > 0, corresponds to a one loop vertex of 
arguments {pf , P2 , P'^ ; , 777,4 , "15}. Recalling Ea. H286|) of ApoendixlEl the whole diagram is therefore 
represented by 



V''' = -^^^ I dCSix-z„Z2) 



x(l — x) 



-e/2 



Jl 



= ^ |2 {Auv - h Coo + Coo dx In A;/ + Inx + ln(l - x 
+ / dCS{x- zi,Z2) ln(l + Q^.;,)|. 

J Bi[Zi,Z2) ) 



In Ea.(|88|) we have a one-loop three-point function to be evaluated up to O (e) (A^^ is defined in Eq.Q). 
The corresponding BT factor is as in Ea. (|74() and, therefore, this function can be computed according 
to the methods of Sect. 4 of or by using the new result of Appendix^ In any case, the behavior for 
6j — > is under control. 



121 



6.3 Wave function renormalization: S, 

To obtain an expression for the derivative of S we follow Eqs.(|7ni)-(I771) with vl = I, vl = and 
^4 = /75, i.e. with the replacement 



Xi{x,y,zi,Z2) 



ii = 4- f^34 zi + (^3 - f^l)y + f^l- 



(89) 



However, being bf = \j = 0, it is more convenient to start from Ea. (|83|) since is 2;2-independent; 
introducing 

■Ax;l = fJ'l - Bi{z) = z'^ - fll4^Z + fll, (90) 

we obtain the same result as in Ea. H88() but with C-functions replaced by i?-functions (one- loop two-point 
functions) of the following type: 



B{o;i} = %;i}o - I %;i}i + O (e^) = ^ dz {1 ; z} B-^~'l\ 



B{0;1}l{x) = / dz{l;z}BY\z) In 
Jo 



1 + 





B,{z) 



Ax (1 - z) 



(91) 
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With their help we write 



^121 
Op 



Bio + B 



10 



dx 



In Aj;-! + lnx + ln(l - x) + En 



}• (92) 



Note that Bii^ is well-behaved for B = and that the simultaneous occurrence of ;B = and of Ax = 
or of z = 1 can be treated according to the results of Appendix^ Furthermore, -Bi,{o;i} are generalized 
one-loop two-point functions that have been described in 



7 The y^^^ diagram 

This topology, corresponding to Fig. Efb), can be written as 



4 1^221 _ ,,2e / jn 



^4 jr. 
TT V 



1 



[1][2][3][4][5] 



(93) 



where we have introduced the following notation: 

[l]=ql+ml, [2] = (gi+pi)2 + mi, [3] = {q^ - q^f + ml 

[4] = (92 + Pi)' + mi [5] = (92 + P? + ml (94) 

and where P = pi +P2 (with all momenta are flowing inward). The related Landau equations are given in 
Appendix^ This diagram, not belonging to the y^^^-class, will represent the first example of application 
of a smoothness algorithms different from the BT one. 

7.1 Evaluation of ^^21 

In evaluating this diagram the first step consists in combining propagators [1] — [3] of Eq. (|94)) with 
Feynman parameters xi,X2, 



.n- . r (3) / / ^.<, / .5, 



n + Q'i)^' 



We have introduced the xi^2-dependent quantities 

Rx = (l- xi) pi - X2 92, QI = xi {ml - ml) + X2 {ql + m| - mf) +ml + {l- xi) p{. 
Integrating over qi in Eq. (|95j) gives 



Vr2F221^ .^-e/2^2erfl + ^) / ^82 



X2(l - X2) 



-l-e/2 



d^q2 



(gi + 2 -92 + ^^1)1+^/2 [4] [5] 



where the new {a;}-dependent parameters are 



P. 



l-xi 2 -p{x{+ xi{p{+ m{- ml) + X2{ml- m{) + m2 
Pi = A pi , = 



1 - X2 ' ' ' ' " ^' X2(l-X2) 

Next we combine the remaining propagators with Feynman parameters = 1,2. It follows 



(95) 



(96) 



, (97) 



(98) 



7r2y221^-^-e/2^2.p/3^_J dS2{xi,X2) X2{1 - X2) 



X / d"92 y2 [123] + (2/1 - y2) [4] + (1 - yi) [5] 



-l-e/2 
3-e/2 



dS2{yi,y2)yi 



e/2 



(99) 
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In the above equation [123] = g| + 2 P^,- • 92 + Af^. After the g2-integration we obtain the following result: 



g(l,M2^ 
M2 

where Q is defined in Eq. (|lfl|) . Since the diagram is ultraviolet finite we set e = and get 

j dS2{xi,X2) j dS2{yi,y2)XG^- (101) 
where the quadratic form Xg given by 

Xa=X2[Ayl + Byl + Cyiy2 + Dyi + Ey2 + Fy (102) 

and its coefficients are 

A —2 2d— 2 2 — — / 2 2\7~,— 2 /2 2 2\ 

A = X2S2V2^ i> = X SiVi, C = X2X [Sp — Si I'l — S2 1^2 j, = (/U4 — ^5 ~ ■^2 ^^2 J ) 

E = X2{xi (Sp + Si - S2 1^2) +X2{s2l'i - Sp) -X2 + ^J.l/x2), F = X^ fi^, (103) 

where Xi = 1 — xi and x = xi — X2- Furthermore, 

nl = si vl xj + xi {fxj - nl-si vl) + X2 - nl) + nl- (104) 

Note that the yi,y2 integral is exactly a generalized Co-function of Appendix El For internal cross- 
checking we use instead Ea. (|264|) . based on a procedure of numerical differentiation, to be discussed in 
Section [mH 

8 The V^^^ diagram 

The y^^^ topology of Fig. El^d) can be written as the following integral: 

where the propagators are 

[1] =ql+ ml [2] = (qi - ga)' + ml [3] = {ql + mi), 

[4] = (92+ Pi)' + mi [5] = + P)2 + mi, [%] = ql + ml (106) 
If 777-3 7^ m-e then y^^^ is the difference of two V^^^ diagrams, 

V'^^^ = 2 ^^^^(-P^; "71, 7772, 7773, 7774, 7775) - 1/^^^(P^; 7771, 7772, 7776, 7774,7775) , (107) 

777g — 7*7,3 

otherwise it is a special case of G^^^, Ea. (|55|) . with ko = ki = 0, and k2 = pi, k^ = P. This choice 
of the rooting of momenta has the advantage of making the zi integration trivial. The corresponding 
quadratic form Xm ^^^^ ~ Xi where we have used zq = y — zi — Z2 — and, moreover, we have 

performed the integration over zi, further renaming Z2-,z^ as ^i,2;2- Since we are only interested in the 
non-trivial case 7773 = 7775 the Landau equations and their solution for V^'^^ are identical to those for 
y^^^,see Appendix El 
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8.1 Evaluation of y^^^ for ma = me 

We proceed with the calculation of V^'^^, always assuming that all external momenta are time-like. 
In this case we use the special set of BT relations which are valid for G^^^ (Section^ when /cq = 0. For 
this diagram we will choose ki = 0,k2 = pi and ks = P. As a consequence the zi integration becomes 
trivial and we obtain 



V 



141 _ 0(2, 



-e/2 



dCS{x;y,zi,Z2) x {1 - x) ' {I - yf/^-' {y - zi) x^/-'{x,y, zi, Z2), (108) 



where G is defined in Eq. (fTn|) and Xm = Xi- The BT-relation that we need is as follows: 



-2-e 



(109) 



«=1,2 



where, similarly to the treatment of the V^^^ case, we have introduced some special combination, where 
the v1 are defined in Eq.Q and in Eq. (|HH) : 



Z\-M = -(1 - i^f - vl) {yl - /i^g) - 2 (i/^ + - ^\) vl, 

Z2;m = (1 - ^1 - I'l) {J^I + nl)+2 {vl - /i^g) vl, 

— I, .2 I ,,2 ,,2^2 ,,2 /I I ,,2 ,,2^ ^,,2 I ,,2 ,,2^ , ,,4 ,,2 , \ ,,2 



yO;Af 



(z^l + /ig - ^4) - ^35 (1 + - 2^2) {^l + /^3 - M4) + ^35 ^1 + 



(110) 



Furthermore, we have Xm = Aj = A(l,i/i,z^2) moreover 6o;a/ = After performing integration by 
parts we introduce secondary quadratic forms: 



Xi;Af = Xm (^2 = zi) = zl- zi + (/ii - fJ'l)y + fil, 

X2-M = Xm (^2 = Q) = vlzl- {vl + /ig - + (^3 - mI) y + /^^. 

with defined in Eq. (|H4j) . The following BT-relations are available for these functions: 



Xi, 



l-e 



1 



if 1 r 

■ Af V 



{y - l)dy + ^{zi- Zii.M)dzi } X,._ 



where additional BT factors and co- factors (Ea. H269|) have been introduced: 



(111) 



:ii2) 



Z 



ll:M 



/^35 



^1; 



-4 A(i,/ii,/ii). 



Zl2;M 
&2;Af 



2vi 



To write our result in a compact form we introduce Zq-m 

■lnXAf(a;,y, ^1,2^2 



4 1/2 
Am, Z: 



A(l/i,/i3,^4). 



(113) 



^+;A/(2;,y,2:i,2;2) 



1-y 



Af) ^3;Af 



and also 



i;Af -^i+l;Af) 



^Af(2;,y,^) = ln(l - y) - ln(3;) - ln(l - x) - \nXi{y,z). 



(114) 



Once more we integrate by parts and, after a Laurent expansion around e = 0, we obtain the following 
expression for the single ultraviolet pole: 



141 

SP ' 
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1 r 3 f 
'^If = 7 — - o I dS2 Zi-uiv) lnXM(y> 



i=0 

Zl 



1 X -V / 12 



l-ZH;„ln/xi + y (iC72(l + Zii;M-2y) lnxM([l,y,0]i) 



where we have used A^/y from Eq.@ and 

2o-Miy) = Am - Zi-M, Zi^2;M{y) = Kiv - Zx-m, Xuiv^^) = Xm(^i ^.y,^)- 

For the finite part the result is as fohows: 



(115) 



(116) 



1 



dC^ (x ; y, {z}) V4;, 



y" dC5 (x ; y, . 



4MHo;m 

Collecting all terms the final answer can be cast into the following form: 
Am 



Z) V3;M - / dC2 V2;M - Vo;A, . (117) 



&0;ii 



1^8 Am IhXmC^j y> 2:1, 2:2) -^"iZx-^M^zx) L+.m{x, y, zi,Z2) 



z. 



V3;jvi = ^ ^r,M{z) L+.M{x,y; [y,z,0]i) + 2Xnt lnxM{x,y; [y,z,0]i 



i=0 

A A/ 



+ |lnXM(y'2) [4(2y - 1) Am - 4Zi;m - 3 Zi;M(y) -Lm(x, y, z) 



ZL 



bo-M 

+ Am (1 - y) lnxM{x,y,y,z)\ + (22; - l - ZH;M)-^^+;M(2;,y ; [y,z,o]i 

+ 3 lnxM(a;,y ; [y,z,0]i 



4=1 



V2;m= fj] ^.;M(y)^M(x,l; [l,y,0]i) + 2 ^ AM(l-y)l lnxM([l,y,0],) 



i=0 



i=l 



+ -r^ Zii;M-Z^+;M(2;,y,0,0) + (2y - 1 - Zh;m)-^^m(2;,1; [l,y,0]i) 



i=l 



VO; 



X lnXM([l>y>0]i) - ^ii;M lnxM(2;,y,0,0) + 2 (1 - y) lnxM([l,y,0]i) 
Am) + 2^-7 + Zii;M In^^s (2 - In^gJ 



i=l 



;ii8) 



The method will fail for bo-M = but also for bo-M 7^ and A(z^f , ;u|, /x|) = and/or X{1, fi^, ix^) = 
which are non-leading Landau singularities representing normal and pseudo thresholds for reduced T/^^^- 
diagrams. To occur, they require either pi or P to be time-like. 

8.2 Evaluation of F^^^: method II 

When bo-M = 0, with bo-M defined in Eq. (|ll()j) . method I as described in the previous section cannot 
be applied. This is not yet a sign that the diagram is singular, since a singularity will appear only if the 
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point of coordinates Zi = Zi fEa. ()82() ') is internal to the integration domain. From general arguments 
presented in Section HTl see in particular Ea. (|53|) . we know that in this case V^^^ ~ ^/bo-M for bo-m — > 0. 
Therefore the correct procedure amounts to applying the BT algorithm of Appendix ^ only once; after 
that we change y — > 1 — y and carry out the y-integration. Some care is needed because a pole at e = 
is hidden in the parametric integration: for this case we use Eg . (12741) of Appendix O As a result we 
obtain an integral representation with the correct asymptotic behavior which requires the introduction 
of generalized one-loop two- and three-point functions. After putting Xm = (1 — ?/) + B{zi, Z2) these 
functions are given in the following list: 



C{0;i}(") = C'{0;i}o - aeC'{0;i}i + O (e^) = y dS2 {I ; zi} B' ^ ''[zi,z2), 
and they correspond to the generalized functions of Appendix IeI Furthermore we define 

B{zi,Z2) 



C{0;1}1 



^■^2 ^7 T In 



1 + 



A, 



(119) 



(120) 



which is well-behaved for B = and where the case B = = will be treated according to Appendix lUl 
Finally we introduce 



^io;l}(^)=^|0;l|0-/?e^io;l|l+0(^ 



'{0;1}0 



{0;1}1 



dz{l; z}B-^-^'{[l,z,0] 



:i2i) 



which are the generalized two-point functions presented in Eq. (10) of 3^ (with a = — 1 — e) and also 



B 



{0;1}L 



dz{l; z}B-\[l,z,0]^) In 



Anr 



(122) 



which are also well behaved for B = 0. Collecting the various terms we obtain: 



V 
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1 



dx 



Vsp: 



.1 



SP;M 



&0;mM4 

where the residue of the ultraviolet pole (A^^ is defined in Eq.®) and the finite part are given by 

2 



(123) 



"00 -^lO 



VsP;M — Cio — Zi-M Coo + Z^. 

i=l 

1^/;*/ = (Coo — Clo) + 2 [-^1;m Cql — Xm Cil + Zi-m Cqi — Am Ch 

Bql - + i*oi - i^ii 



ln(X A^) (Zi;m Coo - Am Cio)\ - - Zr^^ 

i=l 



1 



+ HXA,) {Bio - Bio)] +^{j ln[A, (1 - z,)+B{zi,Z2 

1 r ^ 

-9 / dCiY.Zr^,^ln[A,il-z,) + Bi[l,z,0]i)]] , 



(124) 



and X = X {1 — x) and '+' refers to the subtraction at A^ = 0. 
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8.3 Wave function renormalization: S, 
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The derivative of S^^^ follows from Eg . (17311 with i/^ = 1, i/| = and fJ-i = ^J-^, i.e. with the replacement 
Xi ~^ & of Ea. (|5^ . Therefore, we write 
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Jo 



-e/2 



dxS^-^'ix), S^-'\x)=G{2,NP) I dS2 x{l-x) ' {l-yYn-^z{y-z)Ci 



— 2-e 



(125) 

where Q is defined in Eq. (|lfl|l . The most convenient way of evaluating this integral is to change variable, 
y ^ 1 — y', and to write 



^, = Accy + Bi{z), Bi{z) =hz^ + 2kz + l, \ = hl-k^, Z 



k 
h' 



(126) 



Next we use one BT-iteration, incrementing the power form — 2— e to —1—e, and the results of Apoendixini 
With 



Bn = BnO + ^ Bnl + O {^) = [' dz z" B] 

^ Jo 

B,{z) 



l-e/2 



Bndx) = I dzz^'ByHz) In 

'0 



1 + 



A (1 - z) 



(127) 



we obtain 



A (x) = h{- - Auv) (Z Boo - 2Z Bio + B20) + Z Boi + Bio - Z Bu - S20 + ^ B21) 



h{\zBoL- ZBi^ + \b2i^ + ^ [ dz{z-\z)\ii 

^ ^ -^x Jo ^ 



1 + 



A.(\-z) 
B,{z) 



hi^ZBoo- ZBio + ^B2o) In 



which give our result for this function, with A^^ defined in Eq.@. 



(128) 



9 The V'^^^ diagram 

This topology, depicted in Fig. [Tfe), can be written as 

where we have introduced the following notation for propagators: 

[l]=ql+ml [2] = {qi+Pf+ml [3] = {qi - q2f + ml 

[A]=qj + ml [5] = {q2+pif + ml [6] = [q^ + Pf + ml (130) 

where P = Pi + P2 and all momenta are flowing inward. The corresponding set of Landau equations are 
discussed in Appendix An algorithm to evaluate this diagram is discussed in the following section. 



23 



9.1 Evaluation of y^si 

As a first step in the evaluation of y^^^ we combine propagators [1] — [3] with Feynman parameters 

„4 T/231 



y23l=r(3) ^^^^ j I J dS^ j-^ 

Here we have introduced the x-dependent quantities 



[4] [5] [6] {ql + 2R,.qi + Ql] 
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Ra: = {l-xi)P -X2q2, = xi (m? - m^) + X2 (g^ + m| - m?) + + (1 - xi) P^. (132) 

Integrating over qi gives 

,2e 



2 T/231 



TT^V 



^ni+2 



dSo 



X2{1 - X2) 



-l-e/2 



d^q2 



(g2 + 2P,..g2 + M2)iW2 [4][5][6] 



where new x -dependent parameters are 



P. 



^ ~ p - X p M^- ^^'^ + ml- ml) + X2 {ml - ml) + m^ 



I - X2 - ' - X X2 (1 - X2) 

Next we combine the remaining propagators with Feynman parameters y^, i = 1, . . . , 3. It fohows 



(133) 



(134) 



^2^23i^^if!lrf4 + l 



dS2ixi,X2) X2{1 - X2) 



-l-e/2 



W2" V" ■ 2 

X / d^q2 \y3 [123] + (2/2 - ys) [4] + (yi - 2/2) [5] + (1 - yi) [6] 



dS3{yi,y2,y3)y^ 

-4-e/2 



^/2 



(135) 



In the above equation, [123] = q^ + 'i Px-q2 + M^. After the g2-iiitegration we obtain the following result: 



V 



231 



a(2,M2) 



M4 



-l-e/2 



dSs{{y})yfu^ 



e/2 rr-2-e 
if ; 



dS'2({x}) X2(l-X2) 

where we use Eq.Q and where the quadratic form Uk is given by Uk = y^ H y + 2 y + L or 

Uj, = Ayj + By^ + Cyj + Dyiy2 + Ey^y^ + Fya^s + Gyi + Iy2 + Jy3 + N, 
with coefficients 

^4 = 52 1^2' = ■Sl ^1 1 = Sp, D = Sp — Siuf — S2 ^2, 

-E = X(-Sp + Si l/f - S2 J^l), F = X(-Sp - Si Z^l + S2l^|), G = -S2Z^| + /i5 - /ig, 
I = -Sp + S2l'i + fJ-l- f4^ J = "^X Sp + fil - flj, N = fil, 



(136) 



(137) 



(138) 



where = /i? . It is easily seen that the matrix H of Eq. ()137p is singular and we may change 
variables yi = y[ + X y^, y2 = y2 + X y^ and ys = y'^, with P^ = X P. The transformed quadratic form 
becomes 

U'^ = yiHryr + 2Kiyr + fy3 + L, (139) 

where y* = (2/1,^2)) and the reduced matrix if^ has elements hij with diagonal elements S2 1^2 > ■^i ^1 ^'^d 
off-diagonal elements (sp — si — S2 1^2 Moreover we have 



.2 ,.2 , „ \ , ,.2 ,.2 r _ ,,2 



Sp + S2Z^|+/^l-Ai5> 



f = X {fi^- + Sp) + fi^- fi^, L = Hq. 



(140) 
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Due to the singular nature of the original matrix H, we have been able to confine the xi, X2-dependence 
in the term linear in and there are no terms proportional to 1/3 or yiy3,y2y3- In this case we can use 
the following relation 



1 + 



1 



1 + e 



2/3 



d 

dy3 



i=l 



~2-e 



bK U-'-' 



with a vector Yi 



Hj. ^ Kr- As a result 6^ is xi, j;2-independent, 

bK = fJ-l- (/ill k2 + /122 kj-2 hi2 ki ^2) G 



-1 
12 



(141) 



(142) 



where Hr ^ h, Kr = k and G12 = /in /i22 — /ii2 is the usual Gram determinant. After the transformation 
the innermost integral appearing in Ea. (|136p is written as 



/ dysvl' _ /_ 

Jo JXys JXy-i 



dy2Uj 



-2-e 



where, as usual, X = \ — X and where we can increment the exponent from — 2 — e to 
xi, X2-independent BT factor. The original quadratic form will be denoted by 

Qo;K{yi,y2,y?,) = C/x(yi,y2,y3), 



(143) 
e with a 

(144) 



but after integration by parts the result contains 6 new quadratic forms. One is the original quadric after 
the transformation, 

Qi;K{yi,y2, ys) = U'^iyi, y2, ya), (145) 

while the remaining 5 arise from surface terms. Among them a special role is played by 



Q2;if(yi,y2) = Sp?/i + {-sp + fil - f4)y2 + + sp) + fil- nl 

= H2yl + 2K2y2 + F2 yi + L2. 



2/1+^6 



(146) 



Indeed the quadratic form Q2-k is also incomplete - the matrix of the quadratic part is singular - and 
moreover the xi, j;2-dependent part is confined in the coefficient of yi. Finally we have 



Qi-K{yi,y2) =y* Hiy + 2Kly + Li, y* = (7/1,2/2), 



3, 



,6. 



(147) 



Our strategy will be as follows: the term proportional to Qi-k can be transformed according to Eq. (|141|) 
with 1 + e — > e, i.e. from power — 1 — e to power — e. Furthermore, the form of C/2;a' makes it possible 
to increment once again its power with a BT factor which is xi, X2-independent; this is possible because 
when we have a quadratic in two variables of the form V = hz^ + 2kiz + 2 k2y + / then we may use 



1 



1 



/i+1 ^ 2(/x + l) 



^M+l 



hl-k 



h 



(148) 



The remaining quadratic forms in two variables, from 3 to 6, contain all terms and Ea. H148|) is not active. 
Therefore, the strategy will be to transform all double integrations into the standard form, yi S [0, 1] and 
2/2 € [0, yi] and to use suitable integral representations for the corresponding generalized C-functions. 
The final result will be as follows: 



V 



231 



^ J dS2{xi,X2) I\ 
i=l 



(i) 
K ■ 



(149) 



25 



For the first term, where the BT-algorithm can be applied twice, we expand around e = and transform 
back the integration variables to the standard simplexes, obtaining 



(1) 



1 



1 /" .„ -^(1;A:) 1 



The functions I of Ea. H15U|) are given by 



k=2 



.{1;3) 



-4 In QO. 



1 



(150) 



(151) 



.{1;2) 



X (n - Y2) In Qf.^ + [x{Yi -yiX)- Yi\ In Q\.,, 
+ X{y2X + Y2) In Ql.^ + X (Y^ + X ya) 

+ ^^y^-Y2 + x[{X -2)yi+Y2]}lnQl.,, + {y2 + Y2 + X [(X - 2) 2/2 - 1^2] } In Q(.^ 

+ {y2 - >1 + X [(X - 2) y2 + n] } In Q?.^ (152) 



7(1) 

d K 



Yi + X{yiX-Yi) hiQ'+{Yi-y2-X 



(X-2)2/2+l^i]} InQf.^, 



where Yi = 1 — Yi etc. Note that, in Ea. (|15(J() the '+ '-distribution only applies to the loj 
Integration by parts has introduced new quadratic forms: 



(153) 
>;arithms. 



Q?;K = ^^Qo;x(yi,y2,y3), 

Q\.K = C^Qi,^{l-XyuXy2,yi), 



-CxQi;K(l-Xy2,l-Xyi,y2), 
Q{.^, = CxQi;^(Xyi,Xy2,y2), 0-% 



CxQi;k{^ - Xyi,l - Xyi,y2), 
Ql^^ = CxQi;K{l-Xyi,Xy2,y2), 
Ql.^ = CxQi.^(Xyi,Xyi,y2), 
K - Cx Qi-k{X y2,X y2,yi), 



Cx = X2 {1 — X2)'^ ■ The rearrangement of terms in Ea. (|15Uj) makes evident that no particular 
arises at X2 0. For the second term in Ea. H149p we have from Ea. (|146l) 



Y' 



K2 
H2 



6' 



^2 
H2' 



I 



(2) 



X, 



(2) 



X2 



r(2) 

I — X2 



(154) 
problem 

(155) 



and, correspondingly 



.(2) 



lYi-Y2 
4 6k- 



-I 



1 + 3 / dS: 



X In Ql^ + X In 



- J dCi [(1 -Y^ + X yi) In + (X yi + Y^) In Q|^] } 
where the quadratic forms are 



(156) 



Q%^ = CxQiAy2A-Xyi), 

Q2;K=C^xQi;k(2/1,1-X2/i), 



Q2;K = CxQl;K{y2,Xyi), 

QiK = CxQiAyi^Xyi). 



(157) 
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For the last term in Eg . (11491) we introduce form- factors of the C-family, defined by 



Co;ii;i2= / dS2{l ; x; y} [ax"^ + by"^ + cxy + dx + ey + f - i6 



(158) 



and obtain the following expression: 



r(3) 



X2 



-Yi [xi Co{a) + xCoib)\ - Y2 Coid) xi - 2x' Co(e) + 2 Co{f) -xY^ Co{h) 
Yi Co{i) + Y2 Co(j) + Y2 Coil)] + X2 Yi Co(m) + ^ Cii(a) 

J X2 



_ X 

xi — 

X2 



3 2 

Cii (6) - Cu id) + 2 ^ Cn (e) - 2 Cn (/) + ^ Cn (/i) 
a;2 X2 



- XI Cu{i) -xCuij) + 2x1^ Cuic) + 2 ^ C^ig) -xCu{l)- xi Cuim), (159) 

X2 X2 

where Xj = 1 — Xj and x = xi — X2. Labels from a to m in Ea. ()159() characterize different C form-factors 
and the corresponding expressions for a, - ■ ■ , f are given in Tab. ^ 



label 


a/x2 


b/x2 


C/X2 


d 


e 


f 


a 


xf hu 


x\ h22 


2xf /ll2 


X2(F - 2X2 Xi fi) 


-2X2X1 X2 (T2 


X2 X2 S 


b 


x\ hu 


X^ /l22 


—2x\xh\2 


X2(F - 2X2X1 Pi) 


2 X2 X2 X P2 


X2 x| u; 


c 


x\ hu 


X^ h22 


—2x\xh\2 


-2X2X1X2 Pi 


X2(-F + 2X2 X P2) 


X2 X2 UJ 


d 


X^ h22 


xf /ill 


—2x\xh\2 


X2{F + 2X2 X P2) 


-2 X2 Xi X2 Pi 


X2 X2 W 


e 


xi x^ hu 


Xi x^ /122 


2x\ x^ h\2 


—X2x{F — 2x2 xi ki) 


2 X2 Xi X2 X k2 


x|(F + X2X1 /) 


f 


xi hu 


X\ /l22 


2x1 /112 


—F + 2 X2 xi ki 


2 X2 Xi /C2 


F + X2 Xi / 


g 


x^ hu 


x^ /122 


2x2/li2 


2 X2 X2 X ki 


X2(F + 2X2 X ^2) 


X2 X2 / 


h 


X^ /l22 


x^ /ill 


2x^/112 


X2{F + 2X2 xk2) 


2 X2 X2 X ki 


X2 Xg / 


i 


x\ hu 


X2 /122 


2X1 X2 h\2 


X2{F - 2X2 Xi 0"1 


-2 X2 xi CJ2 


X2 X2 S 


j 


X^ h22 


X2 /ill 


-2X2 X /112 


X2(F + 2X2 X P2) 


-2X2 X2 Pi 


X2 X2 UJ 


1 


X2 hu 


x^ /122 


-2X2 X /112 


-2x2x1 Pi 


X2(-F + 2X2 XP2) 


X2 X2 W 


m 


x\ /122 


xf /ill 


2X1 X2 h\2 


-2 X2 xl (72 


X2(-F - 2X2 Xi 0"l) 


X2 x| E 



Table 1: Parameters for the C- functions arising from Ea. H158|l . 
We have introduced the auxiliary quantities 



h 



11 



2 I, 2 I, ^ / 2 2\ 

-82^2^ /l22 = -Sll^l, hi2 = -{Sp- SiVi - S2l'2) 



ki = ^{pl- pI- S2i^i), k2 = ^{pl- pl-Sp + S2J^i), 1 = pI, 



(160) 
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as well as the combinations 



S = /ill + h22 + 2hi2 + 2ki + 2k2 + l, 

0"! = hii + hi2 + ki, (T2 = h22 + hi2 + k2, 

pi = hn + ki, P2 = hi2 + k2, uj = hn + 2ki + l. (161) 

Finally we have defined 

f = X {fil - + sp) - pI + pI = (162) 

One can see that, for any of the above functions, the corresponding determinant is proportional to G12. 
All C-functions are computed according to the procedure of Appendix El including a preliminary sector 
decomposition in order to avoid numerical instabilities at the end points of the xi —X2 integration. Finally, 
note that the sub-leading BT-factor is given by b'j^ = — 1/4 A(sp, /i^, /ig). An alternative representation 
will be studied in the next section. 



9.2 Evaluation of V^-^^: method II 

There is a second integral representation for this diagram where we start from Eq. (|143|) . exchange the 
order of integration and directly perform the integration. After that the yi — y2 interval is mapped 
into the standard triangle < 7/2 !^ yi 1 and the net result for Ij^ is a combination of 10 Cq functions 
with {x} dependent parameters. We write hij for the elements of Hj- and ki for those of K^, see Ea. H139() : 
the ten quadratic forms in two variables are 

Qi{yi,y2) = aiyl + h yl + q yi y2 + diyi + Cj 2/2 + fu (163) 
and the coefficients are given in Tab. El The result is 

^^ = 71T7 / dS2X2lK, (164) 



= xl Co([l - 2]) - XI X2 Co([3 - 4]) - XI X2 Co([5 - 6]) - Co([7 - 8]) + xi x Co([9 - 10]), (165) 
where Co([i — j]) = Co(i) — Co(i). Furthermore we have 



/ = ^l^^^'^^l A(X1, X2) = Z^^ - X2 (1 - X2) /il + X2 (1 - Xl) (/i| - + Spj, 
X2{l-X2) 

vl = -Sp x\ + Xl {-Sp + pj- H2) + X2 {f4 - lA) + f4- (166) 
As a result our second expression for the V'^^^ diagram is 



^^''^~W I ^'^2(xi,x2) ^J^'^^^ I dS2{yi,y2)lK. (167) 

Note that there is no severe problem when A(xi, X2) because, as seen from Eq. (|165| ) and Tab.jJl the 
differences C(){i) — Co(« + 1) vanish in the same limit. For technical details regarding the evaluation of 
Co we refer once more to Appendix IeI 
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i 


a 


b 


C 


d 


e 


f 


1 


X2 xl /ill 


X2 xf /l22 


2X2 xf /ll2 


-2X2X1 X2 CJl + X2 A 


-2X2X1X2(72 


X2 X2 S 


2 


X2 x\ /ill 


X2 xf /l22 


2X2 xf /il2 


-2X2X1 X2(Ti 


-2x2XiX2(T2 


X2 Xg S 


3 


X2 x\ /ill 


X2 X2 /l22 


2 X2 Xl X2 hi2 


-2X2X1 X2 (Tl + X2 A 


-2 X2 x\ (J2 


X2 X2 S 


4 


X2 xl /ill 


X2 X2 /l22 


2 X2 Xl X2 hi2 


-2X2X1 X2 CTi 


-2X2 x| <72 


X2 Xg S 


5 


X2 x\ /l22 


X2 x\ /ill 


2 X2 Xl X2 hi2 


-2X2 X2 0"2 


-2x2XiX2 (Ji +X2 A 


X2 Xg S 


6 


X2 X2 /l22 


X2 xf /ill 


2 X2 Xl X2 hi2 


-2X2X^ (T2 


-2X2X1X2(71 


X2 X2 S 


7 


X2 /ill 


X2 X^ /l22 


2X2 X^ /ll2 


2 X2 X2 X /ci 


2 X2 X2 X /C2 + X2 A 


X2 X2 / 


8 


X2 X^ /ill 


X2 X^ /l22 


2X2 X^ hl2 


2 X2 X2 X ki 


2 X2 X2 X k2 


X2 X2 / 


9 


X2 xf /ill 


X2 xf /l22 


— 2X2X1 x/li2 


-2x2XiX2/5l +X2 A 


2 X2 X2 X p2 


X2 X2 


10 


X2 x\ /ill 


X2 X^ /l22 


-2X2X1 x/li2 


-2X2 Xl ^2 Pi 


2X2X2XP2 +^2 A 


X2 x\ UJ 



Table 2: CocfRcicnts a ■ ■ ■ f for the 10 quadratic forms in two variables of Ea. H165() . Here S = hu + 2hi2 + 
/i22 + 2A;i + 2k2 + I. Furthermore, ai — hn + hi2 + fci, (J2 = /122 + /112 + /e2, pi = /iii + fci, P2 = /ii2 + k2 and 
= /ill + 2 fci + /. 

9.3 Wave function renormalization: 

The derivative 5^^^ is computed from Ea. ()136() when we put P2 = and 7715 = mg, using Eas. H164|) - 
(|1H7|) with P2 = and ms = mg- This is equivalent to say that all Cq functions in Eq. (|l(i5j) have 
vanishing Gram's determinant. We briefly recall the strategy to evaluate a one-loop vertex when the 
Gram determinant is zero : let us write Hiy + 2 Kj y + Li for any of the quadratic forms of Eq. (|1()3|) . 
Let us introduce also 

hi= lim Gi2Bi, Bi = Li - Kl Hr^ K,, Xi = -AiKi, (168) 

where G12 is the Gram's determinant and Aj is the co-determinant of the element Hi f^i. For each bi ^ 
we obtain 




(169) 



where Xio = 1, Xi3 = 0, and Qi{j j + 1) denotes contractions, i.e. 

Qi(ol) = Q,(l,y), Q,(r2) = Q,(y,y), Q,{23) = Qi{y,0). (170) 

When 6j = in Ea. (|169|) we have to distinguish between two sub-cases (the subscript i is left understood), 

a) e — cd/{2a) 7^ where the integration is trivial and b) e — cd/(2a) = where the vertex reduces to 

a combination of generalized two-point functions 3 showing a possible singularity of the form (a / — 
dV4)-V2. 
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10 The V^^^ diagram 

The 1/222 topology of Fig.[7|^f) is a non-planar one, therefore showing a novel feature with respect to all 
other two-loop vertices, namely there are two lines in common between any of the one-loop subdiagrams. 
The expression is 

where we have introduced the following notation for propagators: 

[1] = ql+ml [2] = (gi - + ml [3] = {qi - 92 + Pi? + ra\, 

[4] = {qi - 92 - P2f + ml [5] = ql + ml [6] = {q^ - Vi? + ml (172) 

The corresponding Landau equations are discussed in Appendix U All the techniques that we have 
adopted so far in evaluating two- loop diagrams will have severe drawbacks when applied to ^^22 
sider, for instance, the standard parametrization: 

10.1 y222. parametrization I 

As usually done, in the first step we combine the propagators [1] — [4] of Ea. H172|) with Feynman 
parameters z/j, i = 1, . . . , 3, obtaining 

,2e 



^4 T/221 
TT V 



(173) 



with y-dependent quantities defined by 



R11 



iy2 - y3)pi - {yi - y2 + yz)p2 - ^2^2, Ql = yi [pl + ml - ml) +y2 (fe -pi? -pl + mj- ml) 



+ y3 (2P ■ q2-pl+ pl + ml- m|) + mj. 
Integrating over qi gives 

where the y-dependent internal and external masses are given by 



~2-e/2 



(174) 



(175) 



D /2/3 1N , ,yi yi-y2. 

Py = { IJPl + ( 1 )P2, 

y2 y2 1 - y2 



y 



1 



yl pI - (y2 - ya) + 2 yi {y2 -yz)P ■ P2 



y2 (1 - y2, 

+ yi {pI -ml + ml) + y2 {pf - pi - ml + m|) - ys [pj - pi + mj - ml) + m\ 



(176) 



Next we combine the remaining propagators with Feynman parameters Xi^i = 1,2. After the 92- 
integration, it follows 



V 



222 



g{2) / dSsiiy}) y2(i-y2) 



-2-e/2 



where Q is defined in Ea. ()lU() and the quadratic form is given by 

Uh = Hh X + 2 KIj X + L„ = Axl + Bxl + Cx\X2 + Dxx + Ex2 + -F, 



(177) 



(178) 
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and where the coefficients are 



A 
D 



B 



p2 



2 I 2 



m. 



6' 



E 



C=-2pi-Py 
= 2pi-Py 



/2 



(179) 



The fact that we have a quadratic form in two variables for the internal integration, in a context where 
there are two external momenta, makes H non-singular, in contrast to the planar topology V^^^. The 
consequence is that when we increase the power of the integrand using only internal variables this will, 
inevitably, result into a denominator which depends on external Feynman parameters. A BT functional 
relation of Appendix ^ applied to the quadratic form produces a coefficient bfj which can be written 
in a compact form as 



rir 



pIp^ 



(180) 



-2 { \iPl 

mi y. L 



4m 



- {ml + ml 



m^ 



(Pi + m. 



y 



o 2\ 2 



2mlpi-Py\'-^ K-P^Ml - Plml)X{-plmlml)]. 



(181) 



Note that the solutions for [pi + Py)"^ of the equation 6/^ = correspond to the anomalous threshold (the 
leading Landau singularity) for the vertex with external momenta pi , Py and internal masses squared 

^ if and only if 



Py, ml and rrig. We have real solutions for {pi + Py 



or 



p\ > (m-s + mef, A(- 



-P?.,M?. 



PLml)>Q. 



(182) 



In principle the proper strategy in dealing with y -dependent BT factors consists in deforming the corre- 
sponding integration contour, a technique introduced in [SH]. In practice the high dimensionality of the 
hyper-contour makes it extremely hard to obtain a distortion with the requested causal properties. An 
alternative parametrization will be discussed in the following section. 



10.2 y^^^: parametrization II 

As it is well-known, any two-loop diagram can be cast into the form of an integral representation 
where the kernel is a generalized sunset diagram. Starting from the definition of y222 gqg_ (ji7Hi72|) we 
combine propagators [1] — [2] with a parameter zi, propagators [3] — [4] with Z2 and [5] — [6] with z^. 
Next we introduce auxiliary variables 



zi {pI +ml) + (1 - zi)ml, 
Z3 {pl+ml) + (1 - 2:3)7715, 



^2 {pI 



ml) 



+ (1 



Z2) {pI 



ml), 



-Zl P2, 



Z2P1 - (l - Z2)P2, K3 



-Z3 Pi- 



rn) 



As a second step we change variables, qi ^ qi — Ki and (72 ^ 92 — and combine the qi and qi — q2 
propagators with a parameter x. After the (71-integration we combine the residual propagators with a 
parameter y and carry out the q2 integration obtaining 



V 



222 



g{2) / dC^y'+'/^l-y) \x{l-x) 



-l-t/2 



Q2 y2 + _ + Q') y + M' 



(184) 
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where Q is defined in Ea. (|lU() and where we have introduced the following quantities 



Q = Ki-K2- K-i, M' = Ri- Ki, 



Mt 



x{l — x) 



(185) 



Since there are no ultraviolet poles we can put e = and increase twice the power of the quadratic form 
in y, obtaining 



222 



/ 



Blx{l-x 



■V 



222 



where the explicit expression for V will not be reported explicitly and where the BT factor is 

B„ = \i-Q\M\Ml). 



(186) 



(187) 



This formulation has the advantage that, for fixed {^i}, we could distort the x-integration in order to 
avoid the zeros of and the study of these zeros follows directly from the analysis of S*^^^ We 
assume that all external momenta are time- like and put = — P'^ Vq, with P = pi + p2 : it follows that 

2 2 2 2 , 2 2/ 2 2 i\ /i 2 2\ /i , 2 2\ o 2 

= ^2 ^1 - ^2 + ^2 ^3(^1 - ^2 - 1) ^1^2 - (1 - 1^1 - 2^2) ^1^3 - i^ + '^l - ^2)^223 - 2 ^1 

- (1 - + ul) Z2 + {1- vl - vl) ^3 - yl. (188) 



^1 4 - Wi + ^5 - /^e) ^3 + M5> 



(189) 



Mi 



vl zl - {yi + /xi - [li) Z3 + fi'i + (1 - x) 



x{l — x) ' ^ 
In terms of scaled quantities we may rewrite the BT factor as 



Z2 - Z2 + 



B^ = P''X[vl,-ul,- 



x(l — X , 



Furthermore we have that 



x-H„, 



.(190) 



(191) 



(192) 



x{l — X) 

where X = x{l — x) and numerical evaluation requires distortion across the zeros of bn- 

10.3 The analytical structure of F^^^ 

Before attempting an evaluation of \/222 jg i^iportant to know more about its analytical structure 
and, therefore, we start once again considering the expression for 6^; zeros of that are real and inside 
the integration region represent apparent singularities and are an obstacle for numerical integration. In 
the following we classify their nature: in this representation we have introduced effective squared masses 
which are given by 



Ml = XiZi;P2,'mi,m2) 



X2iz2;P'^,'m4,m3) 



= X3(^3;Pi,"Z5,m6) 



2 2 I / 2 2 I 2\ I 2 

-P2 ^1 + ("^2 - "T-l + P2) ^1 + 

zl + {ml -ml + P^) Z2 + ml, 

2 2 i I 2 2 , 2\ I 2 



(193) 



32 



Being functions of {z}, their sign is not constant and the following inequalities hold: 



M? 


> 


for 


»o < 




— »n m? mo) < 










or 


»o > 


Af- 


-T>o 777,? ?77,o) ^ D 
/^Z? ''''li ''''2/ — ^1 


zi « *C zi ■<! Zi ^ 








or 


Po < 0, 


Af- 


-?7o. 772?. mo) ^ 0, 


Za z^n^ or 




Mn 


> 


for 


< 0, 


A( 


-P'^,m'i,ml) < 0, 












> 






^2a ^ ^2 ^ ^2o 








or 


p2 < 0, 


A(- 


-P2,mi,ml) >0, 


-^2 < Z2a-> or 


Z2 > Z2b 


Mi 


> 


for 


<o, 


A( 


-pf,ml,ml) < 0, 










or 


>o, 


A(- 


-Pi,"i5,me) > 0, 


^3a < ^3 < ^3^), 








or 


15? <o, 


A(- 


2 2 2\ ^ n 


zz < zsa, or 


ZZ > Zsb, 



(194) 

where P = pi+P2 and Zia b 

the roots of = 0. If , M| > then > 0. Clearly, if aU masses are 
such that their squares are positive and > there are no real solutions for x to the equation hn = 0. 
Otherwise, for < and positive (effective) masses squared, the solutions are given by 



(195) 



The condition > is given by Vh{zijZ2,zz) > with 

Vh{zi,Z2,Z3) =plzf + P^ z\ +pIz^ + 2p2 ■ PziZ2 " 2pi ■P2 Z1Z3 - 2pi ■ PZ2Z3 

- 2plzi - 2p2 ■ PZ2 + 2pi ■P2 Z3+pI, 

with P = Pi + P2- Therefore we have > for 

pI>0 and Gu < 0, 
P2 ^ {p2 > 0) and G12 > 0, zi- < zi < zi+ {zi < z\- or z\ > -21+), 



(196) 



1 r 

zi± = — 

P2 



pI-P2- PZ2+P1-P2Z3± {Z2 - Z3) a/gI^ 



(197) 



Finally, if M| < or < there are again no real solutions. Let us assume that M2 > (0 < X < 1) 



and also M| > and S'^ 



-Q > 0. We have to discuss zeros of bu, i.e. 



(Mj) ^ = (M3 ± Sf , M3 = ^JMl S = (198) 

M?. Only x+ lies between 

(199) 



The minimum, for Mj(a;), occurs at x± = Mi/[M2 =t Mi), where M, 
and 1, corresponding to 

K')min= (^1+^2)^ 

There are three distinct possibilities: 

1. the root (Mj)+ is below the minimum, i.e. Mi + M2 — M3 ^ S, therefore bu can never be zero 



2. Only one root is above the minimum, i.e. {S — M3) < (Mi + M2) < {S + M3) , when there are 
two values of x where 6^^ = 0, 



2{S + Ma)^ 



{S + M3f - Ml + M| ± A^/2 (5 + Ms)^ , Mf,Ml^ 



(200) 
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3. Both roots are above the minimum, i.e. (Ml + M2) < {S- Ms), when we have four values of x 
where bu = 0. The new pair of points is given by 

1 



2 {S - M^f 



{S - M^f - Ml + M| ± M^^2 ^ Ml.Ml^ . (201) 



Next we study the imaginary part of \d.U . Rewriting f7 as a function of x we obtain U{x,y]{z}) = 
Un{x, y ; {z})/x{l — x) with 

U^{x, y ■ {z}) = -{yQ^ + Mi) (1 -y)x^+ [-y^ Q^ + y {Q^ + M^ - M| - M|) + M|] x + y M|. (202) 

There are two roots for C/jv = to be denoted x^^^. One possible way of computing ^^22 ^q^j^^j 
distort the integration hyper-contour. Consider the distortion for the points x^: this possibility ceases 
when x^^fj pinch the real x-axis at a point G [0,1]. Given that Uj^{x,y; {z}) is quadratic in x, say 
C/jv = o,x^ + bx + c, this situation will occur when we simultaneously have xZ = x^ = and 6^ = 4ac. 
The condition for coincidence, i.e. x~ = x^, is 



S = M1 + M2 + M3, 5=V^- (203) 
The remaining two conditions, namely Imx^^H = and x^, = x^ = xl = xl,^, require 

Ms 

= M1 + M2 + M3' 

and it can be easily shown that we have a pinch at y = yth- A possible way out would require distorting 
the {z} hyper-contour in order to avoid the pinch at y = until a true singularity appears. However it 
is very time consuming to build an automatized algorithm that accomplish the distortion in a proper way, 
i.e. that avoids crossing of cuts in the logarithms. All the attempts that we have made do not satisfy our 
bounds on the required CPU time for evaluating a two-loop diagram. Nevertheless we will present the 
main ingredients of the complete analysis, since understanding the analytical structure of the diagram 
has a role of its own, not bound by the method actually used in the numerical evaluation. There are also 
other solutions to the condition for coincidence, among which 



S = Mi-M2 + M3, S = (205) 

Furthermore, if 

Ms 

y^'^ = Mi-M, + M3 ' 
we have again Imx^^j^ = and However this solution corresponds to 

M2 



=^ M2 - Ml ' 

which lies outside [0, 1] . Actually the branch points pinch the real x axis in the interval [0, 1] for 

M3 (M3 -M2 + Ml) - 2 (^MiMa ± y/A 
(M3 - M2 + Ml 



(207) 



yi,2 = 777 . ...2 (208) 



where A = — M1M2 (M3 — M2) (M3 -|- Mi). With considerations completely similar to the ones we are 
going to illustrate in the case S = M3 — Mi -|- M2 , it is possible to show that the values of y in Eq. (|208|) 
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are not included in [0, 1]. Finally, the equation x_ = admits also the solutions S = M3 — Mi — M2 
and S = M3 — Ml + M2. If 5 = M3 — Mi — M2, the two branching points of the logarithm coincide when 

or when 

Ms (M3 - Ml - M2) + 2 (M1M2 ± \/a) 

yth = — A = M1M2 (Ms - Ml) (Ms - M2) . (210) 

(Ms - Ml - M2)^ 

Since by definition S > and S < Ms , the solution of Eq. H2U9() is larger than 1 and so it is of no interest 
for our discussion. The two solutions of Eq. (|2inj) give a pinch in 

„ M2 , , 

XL = Xr} = X = X, = — -— , (211) 

" " + M2 - Ml ' ^ ^ 

which lies outside [0, 1]. 

If S = Ms — Ml + M2 the two branching points of the logarithm coincide when 

M3-M1 + M2' ^ ' 

or when 

Ms (Ms - Ml + M2) - 2 (M1M2 ± Va"^ 
(Ms - Ml + M2)^ 

where A = — MiM2(M3 — Mi) (Ms + M2). In this case the solution of Eq. (|212)) corresponds to a pinch 
outside the interval [0, 1] on the x axis. For what concerns the two solutions of Eq. (|21H|) they are 
complex if Ms > Mi, since in this case the quantity A is negative. Let us consider what happens if 
Ml — M2 < Ms < Ml. We parametrize Ms as Ms = Mi — aM2 where a € [0, 1]. As a consequence, the 
solutions of Eq. ()213|) will read 



yth = 777 , ,'2 (213) 



-a (1 - a) M| - (1 + a) M1M2 ± 2^5 



yth = ^''1 '^2r.V'^'^'^^ "^" ^ ^ = aMiMl [Ml + (1 - a) M2] , (214) 



(a - 1)^M 



so that we can see immediately that the solution with the minus sign in front of the square root in 
Eq. (|214j) is negative. It is also possible to check that also the other solution is always negative; in fact 
the condition 

- Q (1 - a) M| - (1 + a) M1M2 + 2V6 < (215) 

is satisfied when (Mi — aM2)^ > 0, which is certainly true. 

Consider now the distortion for the points xl^. The condition for coincidence, i.e. = x^, is 



S = Mi + M2-Ms, S=^/^. (216) 
The remaining two conditions, namely ImxL__R = and Xl = = x^ = x^, require 

^ 0[O,11. (217) 

" Ml + M2 - Ms ^ ^ ^ ^ ^ 
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The condition for the coincidence of the two x~\. solutions is satisfied also by the choices 

5 = Mi-M2-M3, S = M2- Ml- M:i,and S = -{Mi + M2 + M^). (218) 

Since the quantity 5 is defined positive, the last case in Ea. H218|) . requiring a negative S, does not need 
further discussion. If 5 = Mi — M2 — M3, the two solutions to the equation 6^ = are pinched by the 
branch points of the logarithm, in the interval [0, 1] on the x axis, when y assumes one of the two values 

-M3 (Ml -M2- Ms) - 2M1M2 ± 2^/6 , , w X , . 

yi,2 = ^ , (^ = - M3-M1 M3 + M2 M1M2. 219 

(M3 - Ml + M2Y 

Since S* > 0, then M3 < Mi — M2 and this automatically guarantees that the values of y in Ea. (|219() 
are real; in particular 2/2 j corresponding to the choice of the minus sign in front of the square root in 
Ea. (|219() . is easily seen to be negative. The solution yi is negative (or vanishes) if 

- M3 (Ml - M2 - M3) - 2M1M2 + 2^/5 < 0. (220) 

The condition in Ea. H22U|) is satisfied if M|(Mi — M2 — M3)^ > 0, which is certainly true. We can then 
conclude that, if = Mi — M2 — M3, none of the values of y for which we have a pinch of is included 
in the interval [0, 1], and so it is not necessary to distort the integration contour on the complex y plane. 
Finally we can repeat exactly the same discussion and reach the same conclusions in the case in which 
S = M2 — Ml — M3; the necessary equations can be obtained from the ones used in the S = Mi — M2 — M3 
case by exchanging Mi and M2. 

The general analysis will be as follows. If > and M3 > or < and M3 < there are no 
real solutions for x to hf, = 0, therefore no distortion is needed. Otherwise we have 



< 0, M| > ^ p+ = +{S± Msf, S = M3 = J Ml 



3 ' 



> 0, M| < ^ p- = -{S± Msf, S = ^/Q^, M3 = ^J-Ml (221) 
In order to have real roots for we have to require 

p^-2(Mf +M|)p+ + (Mi2-M|)2 >0, or pi + 2 {Nf + N^) + {Nf - N^f > 0, (222) 



where we have introduced Nf = —Mf since, as we are going to see in more detail later on, the condition 
M^ = p_ implies (Mi it M|)^ < 0. In the first case we obtain 

Mf M| < 0, or Mf M| > and |p+ < (Mi - M2)^ or p+ > (Mi + Mg)^}, (223) 
with Mi = \J\M1\. In the second case we obtain 

iVfiV|<0, or iVfiV|>0 and \^p- < - {Ni - N2f or p- > - (iVi + A^2)^}, (224) 



with Ni = ^^\N1\. For < 0, M| > we have already covered the case MfM| > 0, as long as all the 

results are understood with Mf 2 l-^i 2!- The last case to be considered is > and M| < 0. Let us 
introduce some additional notation; the quantities S and A'3 are defined by 



S = Vp^ > and N3 = J -Ml > 0, (225) 
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and then the dangerous denominator in the BT algorithm will read 

bH = -4:S'^Nl + {M^ + Nl + S^f . (226) 

The apparent singularities in the numerical integration are encountered when bn vanishes, and that 
happens when 

Ml = = - (A^3 ± Sf , (227) 



which is exactly the condition already found in Eq. (|221j) . If we now impose that the two values of x at 
which Ea. ()227() is satisfied coincide, we find that this happens when 



{Ns ±Sf = - (Mf + M|) ± 2^MlMl. (228) 

Assuming M^Ml > 0, the above equation becomes 

(TVs ±Sf = - (Ml ± Mif . (229) 

Clearly this condition can be satisfied only if the square on the r.h.s. is negative, and then if the effective 
masses Mi and M2 are purely imaginary. For the sake of clarity we introduce then the following quantities: 
iVi = yJ-M'( > and N2 = y/-M'^ > 0, so that Ea. d^ can be rewritten as (5 ± N^f = {Ni ± iVa)^- 
It is then easy to see that the two roots of b^ coincide when 5 assumes one of the following values 

5+ = iVi ± iV2 + N3 or SZ = -{Ni ± N2) + iVg, (230) 

and when S is equal to 

Sl = Ni±N2- iVs or S+ = -{Ni± N2) - N3. (231) 

It is now possible to check where the solutions of Ea. ()227() are located under the assumption that S is 
given by one of the Eqs. H23UI231|) . These solutions are found to be located at 

where just the value of x corresponding to the plus sign in the denominator of Eq. (l232|) stays in the 
interval [0, 1] on the x axis. 

The branch points for the logarithm pinch the x axis exactly at the same point; in fact, with U = 
ax"^ + bx + c, we have that 

a = -(y52-iV|)(l-y), b = -y^ + y {S^ - nI + + N^) - nI c = -y iV|, (233) 

and the discriminant of the quadratic equation [/ = vanishes if 

2/1,2 = l^[S^ + ^1 - (^1 + ^2)=^ ± ^/]^] , ysA = ^[S' + ^1 - (^1 - ^^2)' ± Vf^] , (234) 

where we have introduced 

pi = [{Ni + N2f - 52 - Nl] ^-4Nl p2 = [(iVi - N2f + + iV|] ^ - 4 iV| 5^ . (235) 

Out of this four solutions, just the first two correspond to a pinch in the [0, 1] interval on the x axis and 
they pinch the axis exactly at the value in which the zeros of bn coincide fEa. (|232p l. 
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It is now necessary to check, for each one of the coincidence conditions hsted in Eas. (|23U() - (|231|) . if 
the values of y corresponding to a dangerous pinch on the x axis are included in the interval [0, 1]. Let 
us start by taking into account the solutions S^. 

a) li S = N-i — Ni — N2 the quantities yi and 2/2 defined in Eq. H234|) become 

Since 5 > and N3 > this solution is positive, but always larger than 1. 

b) If = A^i — + -^3 the quantities yi and y2 defined in Eg . (12341) become 



1 



= — , ,,,2 [N3 (A^i -N2 + N3) - 2N1N2 ± 2^ , (237) 



and pi = - (iVs - iV2) N1N2 {N3 + Ni). Since 5 > we have that iV > A''2 - A^i- It is immediately seen 
that if > N2 the values of y in (|237|) are complex, since the quantity pi is negative. Let us see what 
happens when N2 — Ni < < N2- It is convenient to parametrize the quantity A^3 as N3 = N2 — aNi 
where a S [0, 1], so that Eq. (|237|) can be rewritten as 



-Q (1 - a) i Vf - (1 + g) A ^iA^2 ± 2^6 
(1 - af 



yi,2= "^^ ,.^»,.v,..2^^v. ^ 5 = aN2Nf[N2 + {l-a)m]. (238) 



The solution y2, that shows a minus sign in front of the square root is automatically negative. In order 
to see what happens in the case of the solution yi, it is necessary to see when the inequality 



a 



(1 - a) Nl - (1 + a) N1N2 + 2^ < (239) 



is satisfied. It is possible to check that Eq. (|239j) is true if N'^ {N3 + A'^i — A'^2)^ ^ 0, which is certainly the 
case. 

c) For the case S = N3 — N1 + N2 we can apply exactly the same reasoning and reach the same conclusions 
that have been found for the case S = N3 — N2 + Ni, no distortion of the integration contour on the x 
complex plane is needed. The relevant formulas for this case can be simply obtained from the formulas 
of the previous paragraph by exchanging A'^i with A'^2- 

d) If 5 = A^i + A''2 + the quantities yi and y2 defined in Ea. ()234() become 
this solution falls into the interval [0, 1]. 

e) If = A'^i — A'^2 — -^3 the quantities yi and 7/2 defined in Eq. (|234j) become 



1 



yi.2 = — 7772 [-^3 [Ni -N2- iVs) - 2N1N2 ± 2^] , (241) 



and pi = — {N3 — Ni) N1N2 {N3 + A'^2)- The solution y2, corresponding to the minus sign in front of the 
square root in Ea. H241|) . is negative, and does not therefore pose any problem. 

The solution yi requires a more careful analysis; since S" > we have that N3 < Ni — N2. The 
expression under square root is then positive and the solution yi is real. It is then necessary to find out 
under which conditions the numerator of the r.h.s. of Ea. H24ip is negative. It is possible to verify that 
the condition 

- iVa (iVi -N2- N3) - 2N1N2 + 2^ < (242) 
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implies iV| (A^3 + A^2 — -^i) > 0, which is always satisfied; we do not have then to distort the integration 
path to avoid the solution yi . 

f) S = —{Ni + N2 + N^) would require a negative S, which cannot occur as S is positive by definition. 

g) If S" = — -^1 — -^3 we can apply the same considerations written for the case S = Ni — N2 — N3 
exchanging everywhere A'^i with 

10.4 Evaluation of ^^22 

Given the intrinsic difficulty in distorting the integration contour, in this section we introduce our 
alternative algorithm to evaluate 1/^22 _ ^yj^j-^ ^j-^g help of parametrization II, introduced in Section fl0.21 
we can write 

3 

y222^_ dC^{x,y,{z])c,{x,y)Q-\ Qh = c,(x, y) Q2({^}) + ^ y) m2({z}), (243) 

1=1 

where the x, y dependent coefficients are 

Cq = xy{l-x){l-y), ci = xy, C2 = ?/(l-x), C3 = x (1 - x) (1 - y), (244) 
and the {z} dependent ones are 

= Xci^i,^2,Z3), (245) 

= XBiplmlml; zi), M| = Xb (-P^ "il, ; ^2), M| = rn^, "^6 ; ^3)- (246) 

Here we have defined 

Xg{p'^,mf,mj ; z) = -p^ z^ + (j? + m| - mf) z + vr?^, Xci{A) = z*' H z + 2 z + L, 

Hij = ki ■ kj ki = p2, k2 = P, k3 = -pi, Ki = -2p2 ■ h, L = pj. (247) 

Furthermore, we introduce 

Q^ = Ql + S, Ql = z^Hz + 2K^z. (248) 
It follows that the original integral can we written as 

r^2 p 

= d^j dC^{x,y,{z})c^\x,y) In Q% \,^^2, 

3 

Q% = c,{x,y) [Ql{{z]) + s\+Y,Cr{x,y)Mh{z]). (249) 

i=l 

Therefore, for y222^ algorithm is based on numerical differentiation despite the poor reputation 
enjoyed by this branch of numerical analysis, at least on the real axis. For previous applications of 
numerical differentiation in this field we refer to In general, the second derivative of a generic 
function f[x) can be written in terms of a so-called (2 A'^ + l)-point formula j34j 

1 +^ 

/^'^(^) = /^ E c''^f{x + nh) + R^{x), c^„ = c^, 5]c„ = 0. (250) 

n=-N n 
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A typical example is given by A'^ = 2 where cq = 



—30 and ci = 16, C2 = 



— 1. For this case 



R2{x) = —h^f^^\x) + 0{h^). 



(251) 



The crucial point in obtaining a decent estimate of our derivative is the choice of the base interval h 
which is also connected with the estimate of the corresponding error on There are two sources of 

error, one related to discretization which prevents h from being large and one connected with rounding- 
off in evaluating / which prevents h from being too small. In our case, if 6^ f represents the round-off 
error on the 2iV + 1 approximation to /^^^ we select the optimal value for h as the one which minimizes 
f /h'^ + Rn and combine errors from round-off and discretization (e.g. from Eq. (|251|) . which requires 
some estimate of /^^^) in quadrature. Once a value for N is selected we obtain an approximation for y222 
which reads as follows: 



2i,2 

n n 



3 

Ql = Y,c^{x,y)Mf{{z]). (252) 

i=l 

As a final technical detail, we note that we always start by performing a sector decomposition to deal 
with common zeros of Cq and which lie at the vertices of the x, y integration square. 

A potential problem of the approach based on a second derivative is that we have to deal with a 
function possessing integrable singularities but in a five-fold domain. As a consequence it is difficult 
to keep the integration (round-off) error very small and h cannot be chosen too small with an obvious 
effect on discretization error. Therefore we have also considered an additional approach based on third 
derivatives. Since 



Cq{x,y)Q] 



1 r 



CqQ^ + Ql In CqQ^ + Ql 



(253) 



we may use a five-point approximation (or higher) [3^ for the third derivative and obtain after straight- 
forward algebra 



V 



222 



222 

3 I 5 



1 



2h? 



dC^{x,y,{z})c '^{x,y) 



z 1 



h^z'/{l + hzY 

/l2zV(l - hzf 



, , 1 — -ih^z"^ ^ 



where the variable z is defined by 



Cq{Ql+pl) + Ql 



(254) 



(255) 



An higher number of points will give similar results. The advantage of Ea. H254() is that, modulo round-off 
errors, h can be taken to be arbitrarily small. The overall advantage of taking a third derivative is that 
the integrand in Eq. H254|) has the form of function x In (function) , something similar to the original BT 
philosophy, i.e. we can write 



V 



222 



1 



2/l2 



+ 



dC^{x,y,{z})Cq (x,y) 
In 



1-2 



4 In 2 + 2 



--1) In(l-C) 



^ + 1) ln(l+C) 



(256) 
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where = hz. The new kernel and its first derivative are both integrable. The discretization error 
will be proportional to the fifth (or higher) derivative; using a seven-point formula with coefficients 
{-1/2, 2, -5/2, 0, 5/2, -2, 1/2} we obtain the following result 

^3?' = / rfC5(x,y,{z})c-i(x,y)(^Li+L2), (257) 

where Lj = ln(l + Xj) and 

Xi = 48C^(l + 5Xi), ^2 = -60 (1 + 5X2), 
5X, = 20C^ + ^ + O (C^) , = f + y + O (C^) , (258) 

for C — > 0, or 

^^1 + ^2= (3- i) ln(l-3C)+ (3+i) ln(l + 3C) +4(^-2) ln(l-2C) 
-4(^ + 2) ln(l + 2C)5(l-i) ln(l-C)+5(l + i) ln(l + C). (259) 

Eg . (12581) gives the estimate of the error for small C (showing finiteness in the limit C — > 0) while Ea. ()259() 
is more appropriate when (" is finite and the argument of some logarithm crosses zero. Analogous formulae 
where the third derivative expression for 1/^22 -g g^yen in terms of a 7-point (or higher) approximation 
with a discretization error proportional to the 7 th (or higher) derivative have been used and will not be 
reported here. 

Note that sector decomposition is always applied to, say, Eq. (|252j) to factorize the common zeros of 
and of \cq (Qq +^'1) + There are two levels of factorization, a simpler one where we only consider the 

X, y pair, and a more complete one where all variables are taken into account. Details of the procedure 
are given in Appendix lUl 

The same technique can be applied to other diagrams, for instance 

f221 = _ J dC4{x,y,{z})cm{x,y)Qa\ (260) 
3 

Qg = cq{x,y) Q\{z}) + J2 Ciix,y) M^{{z}) - Cm{x,y) mU{z}), (261) 

i=l 

where the coefficients are now 

Cq = x{l-x)y{l-y), Cm = (l-2;)y, ci = 1 - x, C2 = x{l-x)y, C3 = x(l-y), 

= XBiplmlml ; zi), M| = Xs(pi,"^i"^5 ; ^2), M| = ml (262) 
Q\{z}) = (zi pi - Z2P2f -2pizi + 2pi ■P2Z2+ pI (263) 
If we introduce = M^q + S with M^q = —Pizf + (pf + m% — m'{)zi we obtain 

= I '^C'4(x,y,{z}) InQ^I,^^., 
3 

Q% = Cq{x, y) Q^{{z]) + Ci{x, y) Mf{{z]) - c^(x, y) 

i=l 

For the first derivative we will typically use a five-point rule with discretization error proportional to the 
fifth derivative. In Tab. 01 we have shown coefficients for numerical differentiation that are not usually 
found in the literature. 



mU{z})+S 



(264) 
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n/ N 


Co 


Cl 


C2 


C3 


C4 




2/3 


-2 


1 








120 •' 


3/7 





-13/8 


1 


-1/8 




X /^^^ 


4/5 


6 


-4 


1 








5/7 





5/2 


-2 


1/2 






6/7 


-20 


15 


-6 


1 




X /^^^ 


7/9 





-7 


7 


-3 


1/2 





Table 3: Coefficients Ci and leading term R in the discretization error for 2 + 1 - point numerical differentiation 
of a function of one real variable, e.g. Ea. (|250|l . The remaining coefficients for /^"^ are C-i = — ct for n odd and 
C-i = Ci for n even. 

11 Numerical results 

In this section we present numerical results for the two-loop three-point scalar functions. All the 
vertices are evaluated using the Korobov-Conroy number theoretic method with a Monte-Carlo error 
estimate as supplied by the subroutine DOIGDF [36|, a multi-dimensional quadrature, general product 
region. 

The estimated error and the time taken will be approximately proportional to the number of points 
in the chosen Korobov set times the number of random samples to be generated. 

There is no ideal format for presenting numerical results of functions of so many variables and we 
have decided to introduce random tables, i.e. for each entry internal masses, external invariants and 
their signs are generated randomly and then the diagram is computed. Alternatively we have considered 
few physically relevant cases, extracted from processes like Z* ff, H* ff, gg etc. Here the input 
parameter set to be used will be 

M^^ = 80.380 GeV, = 91.1875 GeV, mb = 4GeV, mj = 174.3GeV, Mj^ = 150GeV. (265) 

Kinematic configurations are further defined by the Mandelstam invariants of Eq. (jS)) . 

For most of the topologies we have two methods at our disposal and we perform a comparison; 
however, two results will not be shown when they agree in several digits. Results are shown in Tabs.|SHSl 
and, whenever needed, we fix the ultraviolet pole at e = 1; the unit of mass is set to 1 GeV. 

Clearly, for some topology one method is better than the other and the latter is only used for internal 
cross-check. With so many parameters, ranging over a wide interval of values, it is hard to make an 
overall statement on the goodness of the results which depends on many circumstances as vicinity of a 
n-particle cut [37| or of a region where the diagram changes its sign and it is vanishingly small. Note that 
leading Landau singularities (e.g. the anomalous threshold for the one-loop vertex) do not pose a serious 
obstacle to numerical integration since they are rarely inside the physical region jHE] ; n-particle cuts may 
slow down the numerical evaluation whenever the algorithm used requires only absolute convergence. 

As a general and self-evident rule we observe that whenever the imaginary part of a diagram is zero 
then all algorithms have a practically unlimited precision. To give a quantitative description we consider 
the diagram of Fig. El and evaluate it for different values of the bb invariant mass. Results and numerical 
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Figure 6: A two-loop vertex diagrams of the V^^^ topology for H* bb. 

errors are given in Tab. where the first row gives an invariant mass below any cut; the second row 
above the HH-cut; the third one above the tt-cut and the last one above the Htt-cut. The relative 
numerical error is practically zero when below any cut and 0.40%, 1.06%, 0.43% for the real part and 
0.25% , 0.34% , 0.78% for the imaginary part above the cuts. Relative accuracy for the real part tends to 
deteriorate in the region where it changes its sign while it is increasing for the imaginary part which is 
becoming smaller and smaller. 

For the ^^^2 topology we have been able to compare our results with those of [HHI where the configu- 
ration with Pi=P2 = and = 150 GeV has been considered and where the numerical evaluation has 
been performed using conformal mapping and Fade approximants. The comparison is shown in Tab. ^1 
where the errors of ,39_ have been estimated by comparison of the [8/8] and [9/9] Fade approximants. 
Our estimate of the numerical error is done by adding in quadrature discretization and round-off errors. 

12 Conclusions 

Evaluating scalar diagrams at the multi-loop level is only the beginning of a complex scenario where 
we still have to face many (hard) technical problems before being ready to step forward. Our aim, 
in this paper, has been to set up a collection of algebraic-numerical algorithms to deal with arbitrary 
massive two-loop scalar vertices, therefore adding one more stone to the construction started in The 
main result, therefore, has been to assemble relatively simple expressions for scalar two-loop vertices in 
a systematic and coherent manner so that they can be used for practical calculations. 

As explained in the Introduction, we are going to devote a forthcoming paper to the analysis of tensor 
integrals, i.e. of diagrams with a non-trivial spin structure. Although the main emphasis of this paper has 
been on proving the feasibility of a project for computing fully massive diagrams, one should not forget 
that QED/QCD components are integral parts of the evaluation of any realistic observable. Therefore, 
infrared and also collinear configurations should be treatable within the same class of algorithms, or 
within some extension of them. The corresponding proof is rather lengthy and, also for this reason, we 
decided to have a self-consistent presentation in another paper of this series. 

There are six non-trivial topologies for two-loop vertices, as depicted in Fig.El three of them belonging 
to a special class ~ one-loop self-energy insertion - for which we have given a completely general solution 
in terms of the so-called BT algorithm |2j. The remaining three are, somehow, more difficult to deal 
with. 

Here we should mention that the BT- algorithm is a general approach in the evaluation of multi-loop, 
multi-leg Feynman diagrams, but the general solution for the BT polynomial V of Eq. (|2H6|) is not known; 
even more important, we have fixed a guideline in our project which is based on two principles: an 
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algorithm for numerical evaluation of Feynman diagrams is optimal when there is a minimal number of 
terms in the solution, and the singularities of the diagram are not badly overestimated. 

The first principle is adopted in order to avoid, as much a possible, the problem of cancellations among 
the (usually many) terms in the final answer. Trading one difficult original integral in favor of many easier 
ones with better numerical convergence most often will show a clash with the inherent existence of severe 
numerical cancellations. The second principle is dictated by the need of obtaining meaningful results 
even across thresholds, pseudo-thresholds and anomalous thresholds, and not only in asymptotic regions, 
s small or very large. 

Therefore, for some of the diagrams, we have abandoned the BT-algorithm in favor of alternative 
smoothness algorithms which are defined, in general, in Eq.Q. At the same time we mention that the BT 
formalism becomes much less efficient under the same circumstances when other approaches, noticeably 
integration- by-parts formalism j4Uj . reach the limit of manageability under the present technology in 
algebraic manipulation techniques |41j : it is somehow hard to accept that part of our present limitation 
is related to a poor level of technical handling of large systems of linear equations, however, this really 
represents the bottleneck of many famous approaches. In this respect difference equations could still 
introduce some novelty^. 

Understanding these motivations will hopefully explain our preference for algorithms anchored to 
specific classes of diagrams instead of pursuing the search for some universal treatment (the Holy Grail) . 
Therefore, in our study we introduced new algorithms of smoothness and we heavily employed the pro- 
cedure of extended sector decomposition in order to cure possible numerical instabilities generated by 
subtracted integrands, proving the relevance of the method beyond the usual treatment introduced for 
handling dimensionally regularized end-point singularities, of infrared or collinear nature. 

Our formulation provides a systematic and economical way of evaluating a large number of complicated 
Feynman integrals. To a large extent our approach for solving a multi-scale problem is orthogonal to 
integration-by-parts techniques as applied in recent calculations. The reason is that the latter approach 
is naturally tailored and practically unbeatable for a massless calculation, e.g. for pure QCD/QED, or 
in general for a problem with very few scales, and solutions of recurrence relations for the general setup 
at the two-loop level are poorly known. 

Another requirement that we impose in choosing a computational strategy for any given diagram is 
that it should work in any region, independently of the signs of the Mandelstam invariants and even in the 
unphysical region, e.g. below the production thresholds. Admittedly, this is a rather severe constraint; 
for instance there are algorithms, that we have not presented in this paper, which will describe with high 
accuracy the non-planar ^^22 (^j^gram below the two-particle cut but which fail above it. The one that 
we have presented, although suffering from the drawbacks of numerical differentiation, is rather robust 
in all regions. 

There are many papers in the literature dealing with two-loop vertices in one approximation or the 
other, but almost none presents tables of numerical results and, therefore, we had little material for 
comparison. For this reason we have privileged as much as possible some procedure for internal cross- 
check. 

The complexity of the calculation increases with the number of propagators in a diagram and this is 
also reflected by the time needed for evaluating a six-propagator graph with respect to a four-propagator 
one. However, the CPU time requested by one scalar configuration should not be taken as the unit for 
realistic evaluation of physical observables since in the final procedure entire blocks of diagrams will be 
mapped into a single integral. 

As a consequence, we plan to organize any realistic calculation according to building blocks that are, 

^S. Laporta, private communication and |42|. 
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by construction, gauge-parameter independent and will be computed within our approach. For this we 
need to control the gauge-parameter dependence of individual Green's functions; the tool to be employed 
for this purpose is represented by the use of Nielsen's identities 



Acknowledgments 

We gratefully acknowledge discussions and comparisons with Andrei Davydychev, Misha Kalmykov, 
Volodya Smirnov and Oleg Tarasov. We would like to express our gratitude to Ettore Remiddi for 
important discussions on the evaluation of multi-loop Feynman diagrams numerically. The contribution 
of Stefano Actis to the general project of algebraic- numerical evaluation of Feynman diagrams and to 
several discussions concerning this paper is gratefully acknowledged. The work of A. F. was supported 
by the DFG-Forschergruppe Quantenfeldtheorie, Computeralgebra und Monte-Carlo-Simulation" . 



45 



A Bernstein-Tkachov functional relations 



The Bernstein-Tkachov theorem jSJ tehs us that for any finite set of polynomials Vi{x), where x = 
(xi, . . . ,Xjv) is a vector of Feynman parameters, there exists an identity of the following form: 

V ix,d)l[ Vt"-^ [x)=bW Vt (x) . (266) 

i i 

where P is a polynomial of x and di = d/dxi] B and all coefficients of V are polynomials of /ij and of the 
coefficients of Vi{x). 

For a generic quadric we have an explicit solution for the polynomial V which is due to F. V. Tka- 
chov (see also ref. |44j and ref. (IS]): 



G = j dxV^'ix), (267) 

where the integration region is Xj > 0, Xj < 1 and where V{x) is a quadratic form of x, 

V{x)=x^Hx + 2K^x + L. (268) 
The solution to the problem of determining the polynomial V is as follows: 

V = l- (^ + -^)\^^ ^ X'=K'H-\ B = L-K'H-^K. (269) 
2 (// + 1) 

B is the so-called BT factor and the vector X is usually referred to as the BT co-factor. 



B Complex masses 

In our approach, designed for numerical evaluation of Feynman diagrams, complex masses do not pose 
a problem but some care has to be taken, because a complex pole does not lie on the usual physical (first) 
Riemann sheet. Its location is determined by the fact that it should smoothly approach the value for a 
stable particle when the coupling of the theory tends to zero. The complex pole is rewritten in terms of 
real quantities and F^ as pv = — iV^rriy . Consider for simplicity a scalar one-loop two-point 
function that we compute as 

2 s 

Bo(s) = 7-ln7r-ln^- / dx Inx(x), x{x) = s x'^ + (pv2 - Pvi - s) x + pv,. (270) 

e fJ''^ Jo 

For real masses the correct procedure amounts to replace x with x — ^(^i where S ^ In our case, if 
Tyifrivi > ^v2^v2 then Imx > for < x < 1. Taking into account the correct location of the complex 
poles implies the replacement of Inx with Inx — 2i7r 6{—Rex)- 



C A useful integral 

In this Appendix we consider a type of integral that occurs frequently in this paper and that can be 
evaluated by means of hyper geometric functions. Let us define 

1= [ dyy'/^-^{ay + b-i5)-^-'. (271) 
Jo 
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This integral is divergent in the hmit e — > and some care is needed in its evaluation. As a first step we 
write 

I = 2^^b-'-^Ml + e, 1; 1 + 1; -^y). (272) 

It is more convenient to express this result in a different form, by using well-known properties of hyper- 
geometric functions PT] : 

e 1 e a.., r2(l + e/2) faY\~^/2 



' '^)""%^i(l + ^>l + ^2 + i;-4,). (273) 



e + 2\bJ ' ' 2 2 aY' 

The remaining hypergeometric function is only needed at 0(1) and after collecting all terms we obtain 
the following expansion for the integral: 

/ = ^ _ 1 + _L) + o (e) . (274) 

In this way we have isolated the pole at e = 0. 



D Integrable singularities and sector decomposition 



One of the main problems in numerical multidimensional integration is to handle integrable singu- 
larities lying in arbitrary regions of the integration volume. Our experience in dealing with multi-scale 
Feynman integrands is such that extensions of standard techniques IIS] are to be preferred to procedures 
that automatically adapt themselves to the rate of variation of the integrand at each point. 

To give an example suppose that we have to evaluate 



dx 



1 1 
dy — In 







1 + 



ax + y 



a > 0. 



The integral is well defined as it can be seen after performing the x-integration, 



dy 



Li2 



y 



Lig 



a + 1 

y 



(275) 



(276) 



however what we want is a numerical integration. A source of numerical instabilities is connected to the 
region where x ~ y ~ 0, since both numerator and denominator are vanishing small in the argument of 
the logarithm. A nice solution is to adopt a sector decomposition to factorize their common zero. We 
obtain 



dx 



dy 



In 1 + 



1 



+ - In 

X 



1 + 



X 



(277) 



a + y / X V ax + 1- 
A simple numerical exercise shows that one can gain several orders of magnitude improvement in the 
returned error by using Ea. ()277() instead of Ea. (|275( l. This simple example can be easily generalized to 
more complex situations with many variables, although in this case the number of sectors may increase 
considerably. For special values of external parameters a singularity may develop, for instance 



J(a) 
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dy — In 
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1 + 
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In 1 + 



X + ay 
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l + ay 



+ - In 
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X + a 



(278) 
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which, after the sector decomposition shows that a = is a singularity of J. We are going to ihustrate 
the same technique in a more reahstic example: suppose that one wants to compute 



H 



1 



1 _ 

dx I dy — In 
Q Jo X 



1 + 



ax + x{y) 
< y+ < 1. 



x{y) = h{y- y-) {y - y+) - i 5, 



(279) 



5 — > 0+. Suppose also that < y_ < y+ < 1. First we split the y integral into three intervals, 
[0, y+] and 1], then we change variables according to y = y- y', y = (?/+ — ?/_) y' + y-, and 

y = (1 — y+) y' + y+, respectively. In this way all the zeros of numerator /denominator are located at the 
corners of [0, 1]^ and we can apply a sector decomposition to obtain 7 sectors giving the following result: 
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I x J 



(280) 



Hi = y_ In 1 + 
+ Ay In 1 + 
= y-(i - y) In 
+ (1 - y+) In 



a - hy+ y- (y_ (1 - xy) - y+) 
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a - h{Ayyy 
1 + 



+ (1 - 2/4 
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In 



+ Ay In 
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1 + 
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X 



+ Ay In 



a + h{l- y+Yxy"^ + /lAy (1 - y+)y 

X 



1 + 



ax — h (Ay)^ 



ax + /i(l-y+)((l-y+)y + Ay) 



(281) 



with Ay = y+ — y_ > 0. Once again we have been able to gain a much better numerical stability of the 
integrand. Note that the simple examples given in this Appendix are the prototype on which the realistic 
cases of evaluation via numerical differentiation (Section 110. 4j) or via integral representations with Cq 
kernel (Appendix are patterned. 



E The Co(A)-function 

In this Appendix we consider a special family of integrals that often appears in our calculations and 
that can be easily connected to a one-loop C-function. Define 

Co;ii;i2(A; a.../) = j dS2{l;x;y]V-^-^\x,y), (282) 

V{x, y) = ax^ + 6y^ + cxy + dx + ey + f — i S. (283) 
The connection with a scalar one-loop vertex is through the following identification: 

a = -Pl b = -Pf, c = -2 Pi • P2, 

d = M^- M| + P|, e = - M| + + 2 Pi • P2, / = M|. (284) 

For these functions we can use the full set of results of Sect. 4 of [S| or consider new integral representations. 
The new derivation is essentially similar to the one we have for the scalar one-loop vertex |17] but for 
keeping e 7^ 0. Define a to be a solution of 6 + ca + a = 0, and introduce 

A{y) = {c + 2ab)y + d + ea, B{y) = by^ + ey + f. (285) 
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The total result reads as follows: 



Cn = Cn0--\eCni + O{e^), n = 0,11, 12. 



(286) 



First we transform the variable y according to y ^ y -\- ax, so that V{x, y) = A(y) x + B{y), split the 
integral, 



dx / dy 
Jo 



1 pax pa pi p~a pi 

dx dy = dy dx — dy dx, 

J-ax Jo Jy/a Jo J-y/a 



with a = 1 — a and transform again the variables: y = ay' or y = — ay' . Next we use 

1 



dx{Ax + B) 



(Ax + B)- 



-l-Ae 



(287) 



(288) 



and integrate by parts. New functions are introduced 

Ai{y) = A(ay), A2{y) = A{-ay), Bi{y) = B{ay), B2{y) = B{-ay), (289) 

and also 

Qi,2(y) =^i,2(y) + Si,2(y), Q3,4(y) = ^i,2(y)y + ^1,2(2/), Q5,6(x,y) = ^i,2(y) x + 5i,2(y). (290) 
The result is 



Co,n = / dyCo^niy), 
'0 



(291) 



Cii,n = / dyCl^^{y)+ dy dx Cf^^^{x , y) , Ci2,n = a Cu^n + / dyCi2,niy), (292) 

Jo ' Jo Jv Jo 



Subtracted logarithms will be denoted by 

In" Qi,2(y) = In" QiMv) - In" ^1,2(2/), 
etc, and the various components are 
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ln"+i Qi-y In'^+i Q3 



+ 



a 



C 



12,n 



a 
~A. 



a 

^1 

y 



ln"+i Q5 



a 

To 



ln"+i Qe, 



ln"+i Qi - ln"+^ Q3 



ln"+^ Q2 - y ln"+^ Q4 



ln"+^ Q2 - ln"+^ 



(294) 



In the standard analytical approach Cq would be written as a combination of 12 di-logarithms while 
Ci,2 would be reduced to scalar integrals and expressed in terms of Cq and of various Bq functions 
with the usual appearance of inverse powers of Gi2- The approach here is different and aimed to put 
the integrand in a form that is particular convenient for direct numerical evaluation. The coefficients 
a, . . . , / of Eg ■ ^2H'A\} are usually expressed in terms of masses and momenta which, however, may depend 
on additional Feynman parameters. 
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E.l Recovering the anomalous threshold 

In Eg . (12941) we have introduced extra terms in each of the integrals since their total contribution is 
zero. With Ea. H293() we have achieved that the residue of the poles due to Ai^2 = are zero. There are 
two cases where the Ai are nullified: the common one is due to the fact that there is a < yo ^ 1 such 
that Ai{yo) = 0. However, when masses and momenta depend on external Feynman parameters and our 
C-functions are the kernel in the integral representation for a two-loop diagrams A(y) = aiy + ao may 
become zero because ai = oq = 0. Even in this case our representation holds and we may encounter a real 
singularity only when Ai{y) = Bi{y) = 0, Vy. Suppose that we are considering a one-loop Co-function 



with pf 



-m and mi 



ma 



m,m2 = M. Consider now one of the terms in Ea. (|294() . say 



In Qi/Ai; we have a singularity when the zero of Ai, i.e. 



ay 



d + ea 
c + 2ba 



(295) 



is also a zero of which may occur only if s (s — 4m^ + M^) = 0, the anomalous threshold for this 
configuration. The problem is more involved when masses are function of external Feynman parameters; 
let us consider a simple example which, however, shows all the features of more realistic ones. Suppose 
we have to compute 



dz Co{—w? , — m^ , s ; m , z M , m). 



(296) 



We follow the procedure described above for the Co function and consider again one specific term in 
Eq.dlMl), say InQi/^i: 



Bi = s 



1 .1 ^ 
h ^1 

2 I ^2 7\/r2 — 



^1 = a^T (y - yo), 
0' + z' M' -a{y-yo)A, 



1 + 



{c + 2haf Ai 



a {y — yo) s 2am^/3^ + 8m 



6 s + s"^ 



+ a{y 



yo)m^/3^ 



(297) 



where (3"^ = s — Arm? and where A^ = — 2 am? + s and Aiiyo) = 0. If yo £ [0; 1] both numerator and 
denominator in the argument of the logarithm vanish inside the integration domain, i.e. for y = yo and 
z = 0. This is not yet a sign of singularity as it can be seen by using a sector decomposition (as described 
in Sectional) after splitting the y integration interval. For instance, for < y < yo we change variable, 
y = yo (1 — y') and perform a sector decomposition with respect to y, z obtaining 



1 f I 

_ I dy dz In 

a At Jo Jq 



(c+ 26a)2 a A-r 



D 



D = y z'^ M"^ s {ay yo At + (3^ + z"^ y^ M^) 

+ a yo s 2 a rr? + 8 m'^ — 6s m^ + + a y yo m^ + 



(298) 



where only one component has been shown, the one where sector decomposition stops after the first 
iteration. Clearly y = yo, z = does not represent a singularity for I and only At = does; a similar 
analysis applies to all terms in Co- In conclusion, for the numerical evaluation of two-loop diagrams 
which are based on the integral representation of Co just described, a sector decomposition will give a 
much better numerical stability. 

In another example we consider cases where there is no singularity but numerical instabilities may 
occur at the end points of the integration region in the external parameters. Also here some additional 
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work is needed, essentially another sector decomposition is requested. Let us give a simple but realistic 
example: suppose that we have to compute the following integral, 



dxi 



dx2J{xi,X2), J{xi,X2) 



dyi 



yi 



(299) 



where y is a quadratic form in the variables xi and X2 of the kind of Ea. ()283|) . with coefficients 

a = X2{1 - Xi)"^ hii, 6 = 2:2 (1 - Xi)^ /l22, c = 2X2 (1 - Xi)^ /ll2, 



d= (I-X2) 



F-2x2(l-xi)5i 



-2X2(1-Xi)(l-X2)52, f = X2il-X2fS. (300) 



Clearly V = for xi = X2 = 1- A simple sector decomposition gives 



1= dxi / (ix2 X2 J(l — Xi X2, 1 — X2) 

Jo Jo 



In evaluating the internal Co-function in Eg . (12991) we obtain a factorization: 
Ai{y) = X2Ai{y), Bi{y) = xl{l - X2)Bi{y), i -- 
where /i22 + 2 /112 a + hn = and 

Ai^2 = -F ± 2 xi X2 (1 - X2) a (a /122 + /112) xiy - Si - aS2 



1,2, 



(301) 



(302) 



B^ 



S + axiy 



ah22Xiy -2S2 



Bo 



S + axiy 



a /122 xi y + 2 5*2 



(303) 



After that the numerical integration in Ea. H299|) is free from instabilities. 

Our result for the C-functions has the same range of validity of the scalar one-loop vertex, i.e. pi 
and/or p2 and/or P time-like, pi,P2 and P space-like. Indeed the original derivation does not make 
any reference to the actual values of the internal masses, as long as they are real. For complex masses 
and space-like pi,P2 and P the scene changes with respect to the one-loop case since here the effective, 
xi, . . . , Xn-dependent, internal masses do not necessarily have a positive value for the real part of their 
squared values. Once again the presence of an anomalous threshold corresponds to a situation where 
simultaneously a denominator as well as the argument of the corresponding logarithms can be zero. In 
this case sector decomposition should be applied or the methods of ^ should be used. 



F Landau equations for V^'^^ 

Starting with this Section we discuss the Landau equations and their solutions for the six topologies 
that one encounters in two-loop vertices. Their relevance is obvious, as one cannot safely attempt numer- 
ical integration before knowing something about the analytical structure of the diagrams. Furthermore 
these solutions are notoriously hard to derive with standard methods |48j . 

The Landau equations for the V^"^^ topology of Fig. [Tfa) are 

ai {qf + mf ) = 0, 02 {{qi - 92)^ + ml) = 0, 

as {{q2 - P2? + ml) = 0, ((92 - P)^ + mj) = 0, ^""^^^ 

and also 

aigiM + "2(51 - 92)/^ = 0, - a2iqi - q2)ti + a3{q2- P2)^l + ai{q2 - P)f, = 0. (305) 
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The leading Landau singularity occurs for ctj 7^ 0,Vi. We multiply the two equations Ea. (|3U5() by gi^, 
'72^5 P2fj. and respectively. This gives an homogeneous system of eight equations. 



( qi-qi qi- (qi - ^2) \ 

qi ■q2 q2- {qi - g'2) 

qi ■ P2 P2- (qi - q2) 

\qi-P P-{qi-q2)J 



Ql 
02 



/ -qi ■ {qi - 92) qi ■ (92 - P2) qi ■ (92 - P)\ 

-q2 ■ {qi - 92) q2 ■ (92 - P2) 92 • {q2 - P) 

-P2 ■ {qi - q2) P2 ■ {q2 - P2) P2 ■ {q2 - P) 

\-P-{qi-q2) P-{q2-P2) P ■ {q2 - P) J 



If we are looking for a solution where all the Oj are different from zero then the singularity will occur for 



ql 



ql = -fnl-pl + 2q2 ■P2, 



qi- q2 = \ {ml -m\-ml-pl + 2q2-p2) P ■ = \ {P'^ - ml + m\- pl + 2q2 ■ P2). 
After inserting these relations the compatibility between the first two equations requires the condition 



^2 • P2 = 2 - [mi + 7712) + mg). 



(306) 



As a result, it follows that ai = m2lm\ 02. If we use these relations in the next two equations, we obtain 
the following conditions: 



q\ ■P2 



mi 



2[mi + 777-2, 



(pI- {mi + m2)'^ +ml), qi ■ P 



mi 



2(7711 + m2] 



{P"^ - {mi + m2f + ml). (307) 



By inserting these values in the last four equations we obtain the condition for a proper solution; therefore, 
for arbitrary masses, the leading Landau singularity occurs for: 



P' 



1 



2 777^ 



{-pI - {mi + 7772)2 + 777^) {-pj + ml - mj) + 2 ml {pi + pi) ± (Ai X2f'^ 



(308) 



where Ai = X{—pi,ml,m\) and A2 = X{—pl,{mi + 7772) ,7773). 



G Landau equations for V^^^ 

The Landau equations for this topology (see Fig. [7|^c)) are as follows 

Ql {ql + 777f ) = 0, 02 ((^1 - q2? + ml) = 0, 

Q3 (gi + ml) = 0, a4((92+Pi)2 + ml) = 0, a5((^2 + ^')2 + m2) = 0, ^ ' 

and also 

aiqif, + "2(^1 - q2)^l = 0, -a2{qi - 92)^ + "3 ^2/. + "4 (^2 + Pl)^l + "5 (^2 + -P)/. = 0. (310) 

We recall that the leading Landau singularity occurs for / 0,V7. We multiply the two equations 
Ea. ()31U() by qi^, q2^, pifi and respectively. This gives the following homogeneous system of eight 
equations 



[qi 


qi 


qi 


{qi 


-q2)\ 






-qi 


{qi 


- 92) 


qi 


q2 


qi 


{q2 +pi) 


qi 


{q2 + P)\ 




( 


qi 


q2 


q2 


{qi 


- 92) 






-q2 


{qi 


- ^2) 


q2 


q2 


q2 


{q2 +pi) 


q2 


{q2 + P) 




as 


qi 


pi 


pi 


{qi 


- 92) 






-pi 


{qi 


- 92) 


pi 


■ q2 


pi 


{q2+pi) 


pi 


(92 + P) 




a4 


\qi 


■ p 


p 


{qi 


-q2)/ 




V 


-p 


{qi 


- 92) 


p 


q2 


p 


{q2+Pi) 


p 


{q2 + P)j 




\"5/ 



= 0. 
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A proper solution will occur for 



-m 



3' 



q,.q2='^{ml-ml-ml) Pi ■ q2 = \{-pl + ml - ml) P ■ q, = \{-P^ + m^ - m^, 



Compatibility between the first two equations further requires the condition m^ = mi + 7712; as a conse- 
quence, it follows that ai = m2/mi 02- If we use these relations in the next two equations, we further 
obtain 



qi -Pi 



mi 



2{mi + m2 



{-pl + {mi+m2) 



m 



qi-P 



mi 



2{mi + 7712) 



{-P^ + imi + m2)^-ml). (311) 



By putting these values in the last four equations, we finally derive the condition for a proper solution; 
therefore, for general masses, the leading Landau singularity occurs for: 



pI 



1 



2 ml 



-pl + {mi +m2)'^-ml) (-P^ + +^3)^ -m^) + 2 (mi +m2)^ {pl + P^)±{Xi Aa)^/^] (312) 



where Ai = X{— pi, ml, ml) and A2 = A(— P^, mg, m^). 



H Landau equations for V^^^ 

The Landau equations for this topology (see Fig. [7|^b)) are as follows: 

"3 {{qi 



ai {qf + mf) = 0, 
"4 {{q2 + Pi)'^ + ml) 

and also 



0, 



"2 {{qi +Pi)^ + ml) = 0, 
"5 {(q2 + P)^ + ml) = 0, 



q2)^ + ml) 



(313) 



aiqifj. + 02(91 + pi)^, + 03(91 - 92)/^ = 0, 
-03(91 - 92)^* + 04 {q2 + pi)^, + 05 {q2 + P)^, 



(314) 



The leading Landau singularity occurs for / 0,Vi. We multiply the two equations Ea. ()314() by gi^, 
92/i, Pifi and p2^, respectively. This gives an homogeneous system of eight equations 







/ 91 • 91 


91 • 


(91+ Pi) 91 


(91 


-92)\ 








92 • 91 


92 • 


(91+ Pi) 92 


(91 


- 92) 


( 






Pi ■ 91 


Pi ■ 


(91 +Pl) Pi 


(91 


- 92) 








\P2 ■ 91 


P2 ■ 


(91 +Pl) P2 


(91 


-92)y 






-9i • (91 - 


92) 


91 • (92 +P1) 


91 • 


(92 + ^')\ 




-92 • (91 - 


92) 


92 • (92 +?'l) 


92 • 


(92 + P) 




-pi ■ {qi - 


92) 


Pi ■ (92 + P1) 


Pi 


(92 + P) 


V 


-P2 ■ (91 - 


92) 


P2 ■ (92 +Pl) 


P2 


(92 + P)/ 




If all Qj are different from zero, the singularity will occur for 



qf 

91 -Pi 

92 ■P2 



-mf 



mi 



i(m| 



ml - Pi) 
ml +pI - 



ql 

91 • 92 



P2). 



-ml -pl-2q2 -pi, 

-\ {ml -ml + ml+pl + 2q2- pi) 
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Compatibility between the first three and between the S*'^, 6*^^ and 8^^ equation requires the conditions 
q2-pi = -pf + [{p1 + ml- ml) {ml - m| + m|) ± (Ai A2)^/^] 

qi-P2 = -\{P^ - - pi) + 732 [(-pi + "^4 - ml) {ml -ml+ m|) ± (Ai Ag)^/^] (315) 



4 ml 



where Ai = A(m2,m|,m|), A2 = X{—pi,m1,ml) and A3 = A(— ""t-I; "^5)- By inserting these values in 
the remaining equations, we obtain that the compatibility requires: 
1 



pz 



Ami "^4 



— mg + 7714) {pi + m^ — 771,2) (~P2 + "^4 ~ "^5) ~ ^ "^2 "^4 (Pi + P2) (316) 

+ rs {ml -ml + ml) (A2 \z)^'^ + r {pj + mf - ml) (Ai Ag)^/' + s (-p^ + ml- mj) (Ai Az)^/' 
where r,s = ±1. These are the four possibilities for the leading Landau singularity of the diagram. 



I Landau equations for V^^^ 

The Landau equation for this topology (see Fig. [Tfe)) are as follows: 



Qi {qf + mj) = 0, 

04 (g| + 777|) = 0, 



"2 {{qi + P)^ + ml) = 0, as {{qi - ^2)^ + mj) 



"5 ((^2 +Pi) + w,^ 



0, 



"6 ((92 + P)' + m'e 



■0, 
0, 



(317) 



and also 



"igi/i + 02(91 + P)f, + "3 (91 - q2)^l = 0, 

-Q3(gi - ^2)/^ + "4 92/^ + "5 (0^2 + Pl)fM + "6 (^2 + P)fM = 0. (318) 

The leading Landau singularity occurs for ^ 0,Vi. We multiply the two equations Eq. (|318j) by gi^, 
q2^i, Pi^i and Pfj_, respectively. This gives an homogeneous system of eight equations. If all Oj are different 
from zero we may use 



qi 



2 

-771-1 



P = 1 (-P2 + ^2 



771 n 



1l 

q2 



-771 



4; 



■Pi = \ {-Pl+ml 



m^ 



qi ■ q2 



q2 



777,1 



777^ 



. P = 1 (-P2 + ^2 



777r 



Compatibility requires first of all that 



1 



2 77l2 



(^2 _ ^2 _ ^2) (^2 _ ^2 _ ^2) ^ 2 7771 (m^ + 777^) ± (Ai A2)^/2 



(319) 



where Ai = A(777f , 7773, ?77|) and A2 = A(777|, 777^, 777g). By inserting back this result into the system we 
obtain the following condition: 

1 



qi -Pi 



Ami 



[m^ - 7773 -I- 7774) {-Px 



777. 



777.^ 



±(Ai A3)i/2 



(320) 



A proper solution requires yet another relation among the physical parameters. Therefore we have four 
possibilities in searching for the leading Landau singularity of the diagram: 
1 



vl 



2m\ 



[m^ 



777q 



777. 



[mo 



777, 



777; 



\)+2ml {ml + ml) + {X^ X^)'/' 



1 



4777.2 7774 



r 2 2 2\/2 2 2\/2, 2, 2\ a 2 2 r 2 , 2\ 

[m-^ — 7773 — 7774) (7772 ~ ^-3 ~ "^eJ iPl + "^4 + ~ ^ ^3 "^4 ("^5 + "^6/ 



±{ml -ml- ml) (A2 X^)'/' T {ml -ml- m^ (Ai X^)'/' - {pj + ml + m^) (Ai X2)'^' 
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2 ml 



Ami "^4 



{m\ — ml — ml) {m\ — ml — m^) + 2 ml {m\ + mg) — (Ai A2)^^^ 



[mi — m^ — m^) (m2 — m^ — m^) (pi + 171,4 + ^-sJ ~ 4 7713 777.4 ("^5 + "'-eJ 



^{mj -ml- ml) (A2 As)^/^ T (ml -ml- mj) (Ai Ag)^/^ + {pj + ml + ml) (Ai A2 



,1/2 



(321) 



where we have defined A3 = X{—p1,ml,ml). The upper(lower) sign refers to the +(— ) sign in eq. 
EQ. (|m . 



J Landau equations for V^"^"^ 

The Landau equations for this topology (see Fig. [7|^f)) are as follows: 



ai {qf + r7T,f ) = 
"4 {{qi - q2 



P2)^ + ml) 



0, 



"2 ((gi -P2)^ + "7^) = 0, asiiqi 

"5 {qi + ^5) = o> "6 ((92 



q2+Pi)^ + ml) = 0, 



+ 777^) 



0, 



(322) 



and also 



aiqi^ + a2iqi - P2)^l + as^qi -q2+Pi)fj. + 04(91 - 92 - P2)fM = 0, 
-asiqi - 92 + Pi)tM - Oi4 {qi - 92 - P2)fM + as 92^ + "6(92 - Pl)^J. = 0. 



(323) 



The leading Landau singularity occurs for ai ^ 0,V-i. The strategy for determining the solutions is, 
as usual, to multiply the two equations Ea. H323p by qi^, q2^, pif_i and p2^ respectively. This gives an 
homogeneous system of eight equations. However, finding the general solution of these equations is an 
arduous task and even so we do not learn much and, for v'^'^'^, it is more convenient to study occurrence 
of singularities directly in terms of distortion of the integration hyper-contour. However, some physically 
significant case can be discussed. We split our system into two systems of four equations and derive 
CKi and 05; in this way we obtain two homogeneous systems, 5i,2) each containing three equations with 
unknowns 02,03, 04 and 03,04,00. The two conditions for non trivial solutions are giving raise to the 
same quartic equation in pi ■ qi. If we consider a configuration with 



pI 



pI 



-m 



7772 = mQ 



M, 



mi 



777 (i^ 2, 6), 



(324) 



a solution is, for instance, Pi • 91 = ttt,^ + 1/2 (P^ + M^). Liserting back this relation into our systems we 
may solve for 02 , 03 and derive a non-zero 04 under the condition = — M"^ = — 4 777,^ . 
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K Tables of numerical results 



M{bb) [GeV] 


ReV 


ARey 


ImV 


AlmV 


200 


-0.04913523 


0.5 X 10"^ 








350 


-0.1552182 


0.6227 X 10-'^ 


-0.2036244 


0.5153 X 10-3 


400 


0.0552437 


0.5846 X 10-3 


-0.1546848 


0.5263 X 10-3 


500 


0.0751724 


0.3209 X 10-3 


-0.0330644 


0.2564 X 10-3 



Table 4: Numerical results for the topology V^^^ of Fig.Elas a function of the bb invariant mass. Real and imaginary 
parts are in units of lO"*. A is the estimate of the absolute error. 



Sp/M 


51 /Ml 


S2/M2 


mi 


7712 


ms 




Re 1/121 


Imyi2i 


H* ZZ 




+/2M, 


+/M, 


+/M, 








M„ 


-0.473001 ± 2.4 X 10"^ 





+7^/45 


+/M, 


+/M, 


M„ 


M„ 


M, 


M„ 


-0.472906 ± 1.4 x 10"^ 







+/M, 








M, 


M„ 


-0.472845 ± 2.2 x 10-^ 





+/V8M^ 


+/M, 


+/M, 


M„ 


M„ 


M, 


M„ 


-0.472922 ± 2.1 x 10-^ 





+/V20M^ 


+/M, 


+/M, 


M„ 


M„ 


M, 


M„ 


-0.476661 ± 1.7 X 10-^ 









+/M, 


M„ 




M, 


M„ 


-0.492846 ± 2.8 x 10-^ 


-0.041473 ± 3.0 X IQ-^ 


+/VmM^ 


+/M, 


+/M, 


M„ 


M„ 


M, 


M„ 


-0.460821 ± 5.8 x 10-^ 


-0.097880 ± 9.6 x 10"^ 



Table 5: Numerical results for the topology V^'^^ of Fig.lJfa). All masses are in GeV and Mandelstam invariants 
are defined in Eq.®. The ultraviolet pole is e = 1 and the unit of mass is also 1 GeV. The process to which the 
diagram belongs is specified with input parameters given in Ea. (|265l) . 
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Sp/M 


Si /Ml 


S2/M2 


mi 


ni2 


ni3 


m4 


ni5 


Re 


Im ^131 


Random 




+/95 


-/22 


-/87 


4 


18 


43 


3 


57 


0.00580331(7) 
0.00580328(2) 


-0.00013086(9) 
-0.00013082(2) 


-/58 


-/6 


-/67 


24 


18 


88 


72 


93 


0.001206431(2) 
0.001206427(8) 






Z* bb 




+/100 


+/mb 


-/nib 






M, 


nib 


M, 


0.2480028 X 10"^ 





+/oUU 


+/'mb 


-/nib 






M, 


nib 


M, 


0.12228(6) X 10~2 
0.122334(6) X 10"^ 


0.27345(7) X 10"^ 
0.27334(1) X 10^2 


+/800 




-/nib 






M, 


nib 


M, 


-0.2183(9) X 10-3 
-0.21744(6) X 10^3 


0.84447(12) X lO^^ 
0.84378(8) X 10"^ 


Z* -^tt 




+/500 


+/mt 


-/nit 








mt 




0.16720(6) X 10"^ 
0.16717(3) X 10"^ 


0.673497(7) x lO^^ 
0.673475(5) x lO^^ 



Table 6: Numerical results for the topology V^^^ of Fig.[7|[c). All masses are in GeV and Mandelstam invariants 
are defined in Eq.JSJ- First entry is obtained with method I of Section |0] second entry is obtained with method 
II of Section in?^ The ultraviolet pole is e = 1 and the unit of mass is also 1 GeV. 'Random' implies that all entries 
are generated randomly, otherwise the process to which the diagram belongs is specified with input parameters 
given in Ea. (|255|) . 
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Sp/M 


si/Mi 


S2/M2 


mi 


7712 






1715 






Random 




+/76 


+/83 


+/90 


43 


50 


80 


35 


89 


0.3269872(1) x 10"^ 
0.3269872(1) x 10"^ 






+/27 


+/6 


-/44 


54 


28 


94 


51 


32 


0.1327195(1) X 10-6 
0.1327195(1) X 10-6 






+/43 


-/52 


-/77 


10 


5 


21 


5 


9 


-0.136824(8) X 10"^ 
-0.13682(4) X 10-5 


0.600728(2) X 10-^ 
0.60072(4) X 10-5 


1 /CO 

+lol 


+/33 


-/95 


97 


15 


56 


93 





0.2209(2) X IQ-'^ 
0.2208360(8) x lO"'^ 


0.53471(4) X 10-6 
0.5347270(6) x 10-^ 


z* -^tt 




+/500 




+/'mt 




M, 




nib 




0.96881(21) X IQ-^ 
0.96866(4) X 10"^ 


-0.131428(23) X 10"'^ 
-0.131401(3) X 10-^ 


+/800 


+/m 


+/mt 












0.4803(13) X 10-8 
0.4815(1) X 10-8 


-0.5335(12) X 10-8 
-0.5326(1) X 10-8 



Table 7: Numerical results for the topology V^'^^ of Fig. dd). All masses are in GeV and Mandelstam invariants 
are defined in Eq.JHJ. First entry is obtained with method I of Section l8. II second entry is obtained with method 
II of Section IH?^ The ultraviolet pole is e = 1 and the unit of mass is also 1 GeV. 'Random' implies that all entries 
are generated randomly, otherwise the process to which the diagram belongs is specified with input parameters 
given in Ea. (|265(l . 
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Sp/M 


si/Mi 


S2/M2 


mi 


m2 


ms 


7714 


777,5 


Re 1/221 


Imy22i 


Random 




+/10.1 


+/1.51 


+/1.52 


1.1 


2.2 


2.5 


4.4 


5.5 


-0.113675(2) 


-0.004701(1) 


-/85 


+/36 


+/86 


6.5 


4.2 


8.2 


5.6 


7.2 


0.002310(2) 


-0.001084(3) 


Z* bb 




-\-/ y o 11 ih 


+/mb 


+/mb 






M, 


■mt 




-0.511(3) X 10-"^ 





+/VlOOOmi, 




+/mb 






M, 


mt 




-0.5253(3) X 10"'' 





+/V5000mfe 


+/mh 


+/mb 




mt 


M, 


mt 




-0.6122(2) X 10"'' 





Z* -^tt 




+/V5mt 


+/mt 


+/'mt 




ruh 


M, 


mb 




0.2232(2) X 10^3 


-0.1819(3) X 10-3 


+/V8mt 


+/mt 


+/mt 




nib 


M, 


rrib 




0.2067(4) X 10-3 


-0.1470(3) X 10-3 


+/V20mt 


+/mt 


+/'mt 




nib 




mb 




0.1638(2) X 10^3 


-0.9395(15) X 10-^ 



Table 8: Numerical results for the topology y^^^ of Fig. H^b). All masses are in GeV and Mandelstam invariants 
are defined in Eq.®. 'Random' implies that all entries are generated randomly, otherwise the process to which 
the diagram belongs is specified with input parameters given in Ea. H2t)5| . 
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Sp/M 


si/Mi 


S2/M2 


mi 


1712 


7713 


m4 


w-5 


me 


Rey23i 


Im y23i 


Random 




-/60 


+/48 


+/29 


7.5 


4.9 


1.6 


5.9 


3.5 


3 


0.1353(3) X 10"^ 


-0.413(7) X 10-6 


-/39 


+/85 


-/93 


3.3 


0.2 


8.8 


7 


7.7 


7 


-0.983(2) X 10"^ 


-0.6672(1) X 10-6 


-/2 


+/86 


+/43 


4.3 


5.7 


4.6 


8.5 


6.2 


2.8 


0.7469(7) X 10-5 




Z* bb 




+/100 


+/mi, 




rn-t 


TUh 




rrib 




rrib 


0.5021(2) X 10"^ 


-0.1633(2) X 10-^ 


+/200 


+/mb 




■nib 


rrih 


M, 


rUb 


M, 


rrib 


0.984(1) X 10-^ 


0.5450(6) X 10-8 


Z* tt 




-1- /a /S m . 

i-/Vomj 


+/mt 




mt 


nit 


M, 


mt 


M, 




0.1487(8) X W-^ 


0.1357(6) X 10-9 


+/^/20mt 


+/mt 


+/mt 


mt 


nit 


M, 


mt 


M, 




0.169(1) X 10-9 


0.185(1) X 10-9 


H* ^bb 




+/200 


+/mb 


+/mb 


M„ 


M„ 


mt 


mt 






-0.4913523 X 10"^ 





+/350 




+/mb 


M„ 


M„ 


mt 


mt 






-0.1560(12) X 10-^ 


-0.2039(8) X 10-8 


+/400 




+/mi, 


M„ 




mt 


mt 






0.554(7) X lO-'^ 


-0.1544(6) X 10-8 


+/500 




+/mb 


M„ 




mt 


mt 




mt 


0.756(4) X 10-9 


0.329(5) X 10-9 



Table 9: Numerical results for the topology V^^^ of Fig.[7fe). All masses are in GeV and Mandelstam invariants 
are defined in Eq.®. 'Random' implies that all entries are generated randomly, otherwise the process to which 
the diagram belongs is specified with input parameters given in Ea. H2t)5|l . 
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AReF222 


Im 1/222 


Aim 1^222 


4.0 


0.73312 
0.7331 


1.8 X 10"*^ 
1.4 X 10"^ 








4.5 


0.61645 
0.6216 


1.3 X 10-^ 
7.8 X 10^3 


-0.33495 
-0.3402 


1.2 X lO-'^ 
7.1 X 10-3 


5.0 


0.51844 
0.5203 


2.6 X 10"*^ 
4.0 X 10-3 


-0.43100 
-0.4442 


2.7 X lO-'^ 
9.3 X 10-3 


8.0 


0.14555 
0.1455 


6.8 X 10-6 
2.0 X 10-3 


-0.5460 
-0.5491 


4.9 X 10-5 
4.0 X 10-3 


20.0 


-0.2047 
-0.2058 


8.0 X 10-^ 
5.4 X 10-^ 


-0.1876 
-0.1864 


3.8 X 10-"^ 
3.7 X 10-^ 


100.0 m2 


-0.0382 
-0.0385 


3.2 X 10-^ 
1.0 X 10-^ 


0.0152 
0.0162 


3.3 X 10-3 
7.1 X 10-5 


400.0 


-0.0036 
-0.00324 


1.3 X 10-2 
3.6 X 10-6 


0.0051 
0.00507 


8.3 X 10-3 
1.9 X 10-5 



Table 10: Comparison with the numerical results of PI] for the topology \/222 Fig.[7fe) in units of 10-^. The 
common mass is m = 150 GeV and p\ = p\ — Q. First entry is from 3SJ, second entry is our result. A is the 
estimate of the absolute error. 
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Figure 7: The irreducible two- loop vertex diagrams V"^'^ (7^0). External momenta are flowing inward. 
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